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ABSTRACT 



The underlying physics of giant radio halos and mini halos in galaxy clusters is still an open question, which becomes more pressing 
with the growing number of detections. In this paper, we explore the possibility that radio-emitting electrons are generated in hadronic 
cosmic ray (CR) proton interactions with ambient thermal protons of the intra-cluster medium. Our CR model derives from cosmo- 
logical hydrodynamical simulations of cluster formation and additionally accounts for CR transport in the form of CR streaming and 
diffusion. This opens the possibility of changing the radio halo luminosity by more than an order of magnitude on a dynamical time 
scale. We build a mock galaxy cluster catalog from the large MultiDark A'-body ACDM simulation by adopting a phenomenological 
gas density model for each cluster based on X-ray measurements that matches Sunyaev-Zel'dovich and X-ray scaling relations and lu- 
minosity function. Using magnetic field strength estimates from Faraday rotation measure studies, our model successfully reproduces 
the observed surface brightness profiles of giant radio halos (Coma, A2163) as well as radio mini-halos (Perseus, Ophiuchus), while 
obeying upper limits on the gamma-ray emission in these clusters. Because of CR streaming transport and the different scalings of the 
X-ray luminosity, Lx, and the Sunyaev-Zel'dovich flux, Y, with respect to gas density, our model is also able to simultaneously repro- 
duce the observed bimodality of radio-loud and radio-quiet clusters at the same Lx as well as the unimodal distribution of radio-halo 
luminosity versus Y; thereby suggesting a physical solution to this apparent contradiction. For a plausible fraction of 10% radio-loud 
clusters, our model matches the NVSS radio-halo luminosity function. Constructing an analytical radio-halo luminosity function, we 
demonstrate the unique prospects for low-frequency radio surveys (such as the LOFAR Tier 1 survey) to detect ~ 3500 radio halos 
back to redshift two and to probe the underlying physics of radio halos. We will make our cosmologically complete multi-frequency 
mock catalogs for the (non-)thermal cluster emission at different redshifts publicly and freely available on-line through the MultiDark 
database. 

Key words. Cosmology: theory - Cosmology: diffuse radiation - Galaxies: clusters: intracluster medium - Radio continuum: galaxies 
- Gamma rays: galaxies: clusters - Astroparticle physics 



1. Introduction 

The presence of large-scale diffuse radio synchrotron emission 
in clusters of galaxies proves the existence of relativistic elec- 
trons and magnet ic fields permeati ng the intra-cluster medium 
(ICM) (see, e.g, f eretti et al.ll2004l) . Diffuse cluster radio emis- 
sion can be se parated in two classes: radio halos an d radio 
relics (see, e.g., Kempner et al.ll2004t iFerrari et al ]|2008). Radio 
(mini-)halos (RHs) are centered on clusters and characterized 
by a regular and unpolarized morphology, resembling the mor- 
phology of the thermal X-ray emission. On the contrary, radio 
relics typically lie at the cluster outskirts, have an irregular mor- 
phology, and often show a high degree of polarization. While 
relic s seem to directly trace structure formation shocks (see, 
e.g., Ivan Weeren et al.l|2010h . the explanation for the RH phe- 
nomenon is challenging and still an open question. 

Two principal models have been proposed to explain 
RHs. In the "hadronic model" the radio emitting elec- 
trons are produced in hadronic cosmic ray (CR) proton in- 
teractions with protons of the ambient thermal ICM, re- 



only a very modest fraction of (at most) a few 



of CR-to-thermal pressure (see Dennison _ 19801: 
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Vest r^ [T982 : Blasi & Colafra ncescol 119991 iDolag & EnBlinI 
2000;'Miniati et al. 2001b a; Pfrommer & En61in'2003','2004a'!^ 
Blasi etal. 2007; Pfrommer et al. 2008; Pfrommer 20oi 
Kushniretal. 2009; Donnertetal. 2010a b; iKeshet & Loebl 
2010: Keshet.2010.: Enfilin et al.. 201 1). CR protons and heav- 
ier nuclei, like electrons, can be accelerated and injected into the 
ICM by structure formation shocks, active galactic nuclei (AGN) 
and galactic winds. Due to their higher masses with respect to 
electrons, CR protons are accelerated more efficiently to rela- 
tivistic energies and are expected to show a ratio of the spectral 
energy flux of CR protons to electrons above 1 GeV of about 
100, similarly to what is observed in our Galaxy ( Schlickeised 
2002). Additionally, CR protons have a radiative cooling time 
larger than that of the electrons by the square of the mass ratio 
and therefore can accumulate in clusters for cosmological times 
(IVolketal.lll"996h . Indeed, CR electrons suffer more severe en- 
ergy losses via synchrotron and inverse Compton emission at 
particle energies £ > 100 MeV, and Coulomb losses below that 
energy range. 

In the "re-acceleration model", RHs are thought to be 
the result of re-acceleration of electrons through interactions 
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with plasma waves during powerful states o f ICM turbulence, 
as a consequence of a cluster merger (see Schlickeiser et al.' 
1987; Giovannini et al. 1993; Gitti et al. 2002; Brunetti et al. 
2004; Brunetti & Blasi 2005; B runetti & L azarian 2007i, 120111: 
Brunetti et al.a2009.) . This, however, requires a sufficiently long- 



lived CR electron population at energies around 100 MeV 
which has to be continuously maintained by re-acceleration at 
a rate faster than th e cooling processes. We refer the reader to 
lEnBlin et al. I dloTTh for a discussion on the strengths and weak- 
nesses of these two models. 

RHs can be divided in two classes. Giant radio halos are typi- 
cally associated with merging clusters and have large extensions, 
e.g., the Coma radio halo has an extension of about 2 Mpc. Radio 
mini-halos are associated with relaxed clusters that harbor a cool 
core and typically extend over a few hundred kilo-parsecs, e.g., 
the Perseus radio mini-halo has an extension of about 0.2 Mpc. 
The observed morphological similarities with the thermal X-ray 
emission suggests that RHs may be of hadronic origin. In fact, 
cool-core clusters (CCCs) are characterized by high thermal X- 
ray emissivities and ICM densities that are more peaked in com- 
parison to ( mer ging ) non cool-core clusters (NCCCs) (see, e.g., 
ICroston et al.ll20()8h . This distinct difference in the ICM density 
structure of CCCs and NCCCs is reflected in the morphology of 
the two observed classes of RHs. 

The RH luminosity see ms to be correlated to the thermal 
X-ray luminosity (see, e.g., iBrunetti et al.l l2009t lEnBlin et alJ 



•ay 

1201 lb . However, a large fraction of clusters does not exhibit 
significant diffuse synchrotron emission. Galaxy clusters with 
the same thermal X-ray luminosity show an apparent bimodal- 
ity with respect to their radio luminosity. Either they harbor a 
RH or they do not have an y detectable diffuse radio emission 
(IBrunetti et al.ll200 9; Enfili n et al.ll201lh . This suggests the exis- 
tence of a switch-on/switch-off mechanism that is able to change 
the radio luminosity by more than one order of magnitude. While 
such a mechanism could be easily realized in the framework 
of the re-acceleration model (IBrunetti et al.l l2009l) . the classi- 
cal hadronic model predicts the presence of RHs in all clusters. 
The failure to reproduce the observed cluster radio-to-X-ray bi- 
modality was one of the main criticisms against the hadronic 
model. Additionally, the classical hadronic model cannot repro- 
duce some spectral features observed in clusters, such as the total 
spectral curvature claimed in the Coma radio halo or the spectral 
steepening observed at the boundary of some RHs. However, the 
recent repo rt of spectral flatte ning with energy of the radio halo 
in A2256 dvan Weeren et al.l 120121) rises the question whether 
such a variability among different source of the same class is 
inherent to the underlying physics or whether it is a signature of 
uncertainty in the algorithms used to pro cess the data. 

In a recent work, Enfilin et al] (1201 ll) assess these problems 
of the classical hadronic model by analyzing CR transport pro- 
cesses within a cluster. While CR advection with turbulently 
driven flows results in centrally enhanced CR profiles, the prop- 
agation in form of CR streaming and diffusion produces a flat- 
tening of CR profiles. Hence, different CR transport phenomena 
may also account for the observed bimodality of the radio lumi- 
nosity in the hadronic model and may also explain the spectral 
features observed in some clusters (^filin et al. 201 1 ). Note that 
these CR transport pheno mena were no t considered in earlier an- 
alytical works (see, e.g.. iPfrommer & Enfilia .2004a) as well as 
in previous hydrodynamic simulations (see, e.g ., Mini ati et al.l 
l2001aHPfrommer et al."2008';'Pinzke & Pfrommer.i20 10). 

More recently, Basu (2012) presented the first scaling rela- 
tions between RH luminosity and Sunyaev-Zel'dovich (SZ) flux 
measurements, using the Planck cluster catalogue. While the 



correlation agrees with previous scaling measurements based on 
X-ray data, there is no indication for a bimodal cluster popula- 
tion dividing clusters into radio-loud and radio-quiet objects at 
fixed SZ flux. While the SZ flux correlates tightly with cluster 
mass, the X-ray luminosity, Lx, exhibits a larger scatter. The 
CCCs predominantly populate the high-Lx tail (at any clus- 
ter temperature) and make up approximately half of the radio- 
quiet objects (Enfilin et al. 201 1). This suggests that the switch- 
on/switch-off mechanism may not operate at fixed Lx but also 
causes an evolution of that quantity. As the cluster relaxes af- 
ter a merger, it cools and forms a denser core. Simultaneously, 
Lx is expected to increase which may simultaneously decrease 
the radio luminosity owing to the decaying turbulence that is re- 
sponsible in maintaining the radio emission in either model (that 
accounts for microscopic CR transport). 

An observational test that is able to disentangle between the 
hadronic and re-acceleration models is the gamma-ray emission 
resulting from neutral pion decays, a secondary product of the 
hadronic CR interaction with protons of the ICM, which is not 
predicted by the re-acceleration model. Such observational ef- 
forts have been undertaken in the last few ye ars (for space-based 



cluster observations in th e GeV-band, s ee Reimer et al.l 12003 : 
Ackermann et al."2010a"b'; 'Zimmer et al."2012'; 'Han et al.''2012 ; 
for ground-based observations in the energy band > 100 GeV, 



see Perkins et al. 


|2006; Perkins 2008; Aharonian et al.||2009ap]; 


Domainko et al. 


2009; ,Galante et al. 2009; Kiuchi et al.1 20091 


Acciari et al.l2009l;lAleksic et al.l2010. 2012b wiflioutbeins able 



to detect cluster gamma-ray emissiorfl Gamma-ray limits en- 
able us to put significant constraints on the relative CR pressure 
contributions. The long observational campaign of the Perseus 
galaxy cluster performed by the Major Atmospheric Gamma- 
ray Cherenkov (MAGIC) telescopes constrains the cluster av- 
era ge CR-to-t hermal pressure to be less than a few percent 
( Aleksic et al. ' 2010, 2012). Comparing the co rresponding up- 
per lim its to predictions of simulations by Pin zke & Pfrommen 
constrains the maximum CR acceleration efficiency at 
structure formation shocks to be < 50%. Alternatively, this 
may indeed suggest the presence of non-negligible CR trans- 
port ^grocessesjnto the outer cluster regions as suggested by 
En fiUn et al.l ( 1201 1() . Constraints at the same level are confirmed 
by the FermZ-Large Area Telescope (LAT) data on the Coma 
cluster (Pinzke et al. 201 1; Zimmer et al. 2012; Han et al. 2012!). 
To summarize, gamma-ray observations are not yet sensitive 
enough to severely challenge hadronic models of RHs. 

An important step towards understanding the generating 
mechanism of RHs could come from detailed RH population 
analyses. To date, we know of 53 clusters that harbor RHs (see 
iFeretti et al.l2012l for an almost up-to-date list). Only two X-ray 
flux-limited studies have been conducted that assess the impor- 
t ant q uestion of the RH f requency in clusters (Giovannini et aT] 
ll999[IVenturi et al.ll2008h . Since the number of RHs in such X- 
ray flux-limited samples is small with typically a few RHs, the 
conclusions on the underlying physical mechanisms of RHs are 
not very robust. Fortunately, this is going to change thanks to the 
next-generation of low-frequency radio observatories such as the 
Low Frequency Array (LOFAR) which officially started opera- 



' Recentlv.'Han et al] i2012h claimed an evidence for diffuse gamma- 
ray emission in the Fermi satellite data of the Virgo cluster which how- 
ever needs to be carefully scrutinized by varying the uncertain fore- 
ground modeling (for a preliminary assessment of the Fermi collabora- 
tion, see Elliot Bloom's talk at the UCLA Dark Matter 2012 conference, 
https ://hepconf .physics.ucla.edu/dm 1 2/) . 
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tions in 20 100 In fact, a deep cluster survey is part of the LOFAR 
science key projects and expected to provide a large number of 
radio-emitting gal ax y clusters up to redshift z ~ I (see, e.g., 
ICassano et al.ll2010l:lRottgering et alJl2012 ). This will hopefully 
permit to clearly determine the RH phenomenology with re- 
spect to galaxy cluster properties such as the fractions of radio- 
loud/quiet, non cool-core/cool-core, and non-merging/merging 
clusters, and to explore the role of different parameters like the 
magnetic field, the CR acceleration efliciency, and CR transport 
properties. 

In this work, we use a sample of galaxy clusters up to z = 1, 
obtained fro m the MultiDark A^-body cosmological simulation 
dPrada et al.ll2012) . We construct a phenomenological model to 
assign a gas density profile to each halo in the MultiDark simu- 
lation that is based on observed X-ray cluster properties. We as- 
sume that the RHs are generated by secondaries of hadronic CR 
interactions with the ICM. In adopting the hadronic scenario, we 
construct a hybrid CR model and merge the resu lt of hydrody- 
namic simulation s by Pinzke & Pfrommerl (|20IO|) and fiie ana- 



Table 1. Number of halos in the final cluster sample 



lytical model by EnBIin et al 



(1201 

inclusion of the latter is particularly important as CR transport 
processes have a dramatic impact on the emission profile and 
normafization. We will show that this approach enables us not 
only to reproduce the radio surface brightness profiles but also 
the radio-to-X-ray and radio-to-SZ scaling relations. Hence, this 
explicitly demonstrates that the hadronic model is an attractive 
model that matches all current observational constraints and pre- 
dicts the luminosity function for the LOFAR cluster survey. 

In Section |2] we explain the methodology of assigning gas 
density profiles to dark matter (DM) halos and demonstrate that 
the model reproduces the observed SZ and X-ray scaling re- 
lations and luminosity function. We then construct the model 
for the CR distribution in clusters. In Section |3] we apply our 
model to observed surface brightness profiles of individual RHs 
and explore the allowable parameter space for CRs and mag- 
netic fields. In Section m we compare the modeled radio-to-X- 
ray and radio-to-SZ scaling relations to current observations and 
show how they vary for diff'erent choices of our CR and mag- 
netic parametrizations. In Section|5] we show the radio luminos- 
ity functions, compare them to current observational constraints 
and provide predictions for the LOFAR cluster survey. Finally, 
in Section |6] we present our conclusions. In this work, the clus- 
ter mass Ma and radius Rt^ are defined with respect to a density 
that is A = 200 or A = 500 times the critical density of the 
Universe. We adopt density parameters of - 0.3, Qa - 0.7 
and today's Hubble constant of H^) - 100 x /170 km s"' Mpc"' 
where h-jd - 0.7. 



2. Methodology 

The two basic ingredients of this work are an A^-body, DM- 
only cosmological simulation that we use in constructing our 
complete cluster sample, and the (non-)thermal emission model. 
Here, we use the MultiDark simulation that will be described in 
Section im along with our final cluster sample. As we will see, 
for any given cluster, the necessary ingredients for the emission 
model are mass, temperature, and gas density distribution. We 
assign to each cluster in that simulation a gas density profile that 
is phenomenologically constructed from state-of-art X-ray ob- 
servations, as will be shown in Section IZ21 We will show that 
this approach can successfully reproduce the known X-ray clus- 
ter characteristics, such as the X-ray luminosity function (XLF), 



redshift z 


number of halos 


0.0 


13763 


0.2 


12398 


0.3 


10783 


0.4 


7789 


0.61 


5187 


0.78 


3372 


LO 


1803 



Note. We show the number of halos in our MultiDark snapshots at red- 
shift z for M200 > 1 X lO'"* r' Mo a 1.4 X 10" /i^^ Mq. 



the luminosity-mass relation, Lx - M, and the Yx - M rela- 
tion, where Yx = M ^^^k^T with an X -ray-derived gas mass Mc,as 
and temperature T (iKravtsov et al.ll2006.) . Moreover, we com- 
pare our Fsz - M relation with SZ-derived measurements. Only 
if our model matches available cluster data on the gas properties, 
we can start to explore different parametrizations of CR physics 
and its implications for the radio and gamma-ray emission, as 
will be explained in Section[ 



2.1. MultiDark simulation and final cluster sample 

The MultiDa rk simulation used in this work is described 
in detail in Pra da et al ] (1201211 It is an A^-body cosmo- 
logical simulat ion done with th e Adaptive-Refinement-Tree 
(ART) code ( Kravtsov et al.lll997h of 2048^ particles within a 
(1000 Mpc /i-i)3 cube. The latest WMAP5 and WMAP7 cos- 
mological parameters were used. This simulation is particularly 
well suited for our purpose because of its large number of clus- 
ters. 

We use the MultiDark halo catalog from its on-line 
databa se^! constructed with the Bound Density Maxima algo- 
rithm (iKlvpin & Holtzmanlll997h . We will mainly use M^qo and 
7?5oo for comparison with existing observatio nal wo rks. We use 
the technique described in Hu & Kravtsovl (1200 3l) to convert 
M200 and R20Q provided by the MultiDark halo catalog to M500 
and /?5oo- In creating our cluster sample we only select distinct 
halos, i.e., those halos that are not sub-halos of any other halo, 
which by definition are not galaxy clusters. 

Additionally, we assume that the main emission mechanism 
in the ICM is thermal bremsstrahlung, whic h is true o nly above a 
particle energy of approximately 2.6 keV (Sarazin 1988). Below 
this energy, there could be other important contributions to the 
emission, e.g., from atomic lines. Therefore, we impose a mass 
cutofM2oo> IxlO'"*/!"' Mo ^ 1 .4x10''' h'l Mp w hich ensures 
k^T > 2.6 keV, assuming the M500 - relation of iMantz et al.l 
(l20T0ah . 

The LOFAR radio observatory is expected to detect RHs up 
to redshift z ~ I. Thus, we make use of different simulation 
snapshots up to z = 1 . In Table [1] we show the total cluster 
number in our final cluster sample at diff'erent redshifts. 



2.2. Gas density modeling 

We decided to use a phenomenological approach and to con- 
struct the gas density profiles directly from X-ray observations. 



www.lofar.org 



3 www.multidark.org. 
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Fig. 1. Electron density profiles of the 31 clusters in the 
REXCESS sample. Grey and black lines represent NCCCs and 
CCCs, respectively. The blue and red lines represent our GNFW 
mean profile for the NCCCs and CCCs, respectively. 



A suitable X-ray sample that provides the needed information 
is the Representative XMM-Newton Cluster Structure Sur vey 
(REXCESS) sample (ICroston et alJ 120081: iPratt et al.l l2009l) . It 
is a sample of 31 galaxy clusters of different dynamical states 
at redshifts 0.06 < z < 0.18 with detailed information on 
the de-projected electron density profiles (Croston et al. 2008). 
In Figure [1] we show the 3 1 electron density profiles of the 
REXCESS sample color-coded by CCCs and NCCCs. 

In order to obtain an electron density profile that we will 
attach to our simulated clusters, we use a generalized Navarro- 
Frank- White (GNFW) profile. 



no 



where x - R/Rc and R^ is the cluster core radius. To reduce 
the dimensionality of our fit, we fix representative values of 
Rc - 0.27?5oo, a = 1 and S - 2.5. We fit the radial density 
profiles of the REXCESS sample in log-log space, separating 
them in the two categories of NCCCs and CCCs as shown in 
Figure [1] We obtain no.NCCC = 1.02 x 10"^ /z^q^ cm"^, no,ccc = 

8.32 X 10-^ h\'f^^ cm-3,/3Nccc = -0.093, and/3ccc = 0.592. The 
resulting fits are shown in blue and red for the NCCC and CCC 
population, respectively. 

The next step is to introduce a mass-scaling in order to apply 
our GNFW profiles to all our clusters. We adopt the gas mass 
fraction-mass scahng, /gas,500 - A/500 of ISun et aTl (l2009l) (and 
adopt their Equation (8)). We can express /gas,500 in the following 
way: 



/gas,. 



gas,500 



Mgas,500 Jo PgasdV 



M500 



M500 



(2) 



with pgisiR) = Pgas = ngifip/iXuXs) where nip is the proton 
mass, Xii - 0.76 is the primordial hydrogen mass fraction and 
Xe = 1.157 the ratio of elec tron-to-hydrog en number densi- 
ties in the fully ionized ICM (ISarazinlll988l) . For each cluster 
i of our sample, we then define a mass-scaled gas profile as 
Pgas,/ = C/Pgas with: 

Ci = (0.0656 ± (0.0064^0) /z7o'-^ 



M500J 



1.04X lOi^/j-jjMoj 



0.135±(0.030j2) 



X™' PgasdV 



(3) 



where gi and g2 are random Gaussian number which we use in 
order to simulate the natural scatter of the gas profiles|3 

Hence, for each cluster in our sample we obtain a gas den- 
sity profile Pgas./ that obeys the observed /^flj,5oo - A/500 relation 
and is uniquely determined by its DM mass Msoo,,- and by the 
property of being a NCCC or CCC. We assign the latter prop- 
erty to every halo depending on its merging history. In particu- 
lar, we make use of the offset parameter Xoff computed for the 
MultiDark halo catalog. This is defined as the distance from the 
halo center to the center of mass in units of the virial radius. 
This parameter assesses the dynamical state of the cluster and 
whether the halo experienced a recent merger or not. Current 
observations reveal a ratio of NCCCs and CCCs of about 50% 
(see, e.g., Chen et al. 2007; S anderson et al.l2()09l) . Since there is 
a correlation between merging clusters and NCCCs, we use the 
median of the Xq^ distribution to separate our sample into CCCs 
and NCCCs (with NCCCs defined to be those halos with the 
larger dynamical offsets). Clearly, this is an over-simplification, 
and future X-ray surveys will have to determine this property 
also as a function of redshift. 

We also account for redshift evolution of the gas profiles. 
While our NCCC and CCC gas profiles as derived from the 
REXCESS cluster sample are merely used to define a profile 
shape, the normalization of the gas profiles is s et by the observa- 
tional /pas..500 - M^Qo rel ation jSun et al.ll2009l) . The 43 clusters 
used in fSun et al.1 (l2009l) have redshifts 0.012 < z < 0.12 with a 
median of z ~ 0.04. Thus, our phenomenological gas profile is 
representative of the cluster population at z » 0. To extend this 
profile to high-z, we include a self-similar scaling of the gas den- 
sity as pgas(z) = E(z)^pg^s(z - 0), where E(z)^ - 0.^(1 +z)^+Oa. 



(1) 2.3. X-ray and SZ scaling relations 



In order to check whether our phenomenologically derived gas 
profiles reproduce the observations, we calculate the b olomet- 
ric X-ray thermal bremsstrahlung luminosity Lboi as in SarazinI 
( ISSSfl and compare our sample result with the observed Lboi - 
M500 relation and XLF0 

To assign a temperature to our model clusters (that is needed 
for calc ulating Lbpi and Yx), we adopt the T - M500 relation by 
iMantzet all (1201 Oah . 



logio 



^B^ci 

keV 



^A-\-B logi, 



( E(z)M: 



500 



lOis/z-g'Mo 



(4) 



The values 0.0064 and 0.03 quoted in Equation l[3j do not represent 
the proper scatter of the /gas.soo - ^500 relation but reflect the parameter 
errors and we rescaled the numerical values to a Hubble constant of hjo 
used in this work. 

^ We check our procedure by fitting each of the 31 REXCESS clus- 
ters with Equation (TJ and calculating Lboi with the measured gas tem- 
perature of each cluster. As a result, we fall short of the observed lu- 
minosity by a mean (median) of about 21% (20%). This is reasonable 
considering that we do not permit the parameters R^., a and S to vary be- 
tween difi'erent objects. Additionally, we neglect atomic line emission 
which may give a noticeable contribution, in particularly for low-mass 
clusters and in the cluster outskirts of larger systems. 

* The mean (median) difi'erence at z = between L(,oi within ^^200 or 
within Rsoo is ~ 5% (~ 7%). While Lboi refers to the quantity calculated 
within Rsoo, we note that the XLF for luminosities calculated within 
R200 will be barely changed. 
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Fig. 2. X-ray and SZ scaling relations. Grey triangles show the MultiDark sample (limited to the mass range covered by observa- 
tions), the black line is the corresponding scaling relation, and the blue line is the observational result. The black crosses represent 
the median values of the quantity in question for a given mass bin (indicated by horizontal error bars), and the vertical error bars 
represent the standard deviation with in a bin. Left . We c ompare the bolometric X-ray luminosity-to-mass relation, Lboi - M500, at 
Z = 0.2 to the observational sample bv lMantz et al. (201Qa) with a median of z ~ 0.2. Center. Yx-Msm scaling relation of our model 
in comparison to the observatio nal sample by Mantz et al. (20 10a) . Right. Fsz - A/500 scaling relation at z = 0.16 in comparison to 
the observational sample by the iPlanck Collabo ration et al. ( 201 ll) with a median redshift of about 0.15. The bottom panels show 
the relative difference to the observational scaling relations. 



where A - 0.91, B - 0.46, and T'd is the c l uster temperature not 
centrally excised (see iMantz et af]l2010ah . iMantz et aP (1201 Oah 
report a scatter of cTyx = 0.060 which we apply to our sample 
using Gaussian deviates. 

In the left panel of Figure |2] we show ho w our model 
Lboi - M 500 relation compares with observations bv lMantz et al.l 
(1201 Oah {all data, see their Table 7). Their sample is com- 
posed of 238 clusters at 0.02 < z < 0.46 with a median 
of z ~ 0.2 and self-consistently takes into account all selec- 
tion effects, covari ances, systemat ic uncertainties and the cluster 
mass fu nction (.Mantz et al. 201 Oh) . For this reason, we com- 



pare the iMantz et al. ( 2010ah date to our model at z = 0.2, 



and limit the comparison to the mass range covered by the ob- 
servations. Overall, there is reassuring agreement between our 
phenomenological model and the data, which probe our model 
most closely on scales around the cluster core radii (which is 
where the contribution to Lx per logarithmic interval in radius, 
dLx/dlogr oc r'ni^Jir) y/k^T, attains its maximum). In Table|2] 
we show our model Lboi - ^^500 scaling relation and its scatter 
for different redshifts. We find that the scatter of our samples at 
all redshifts are Gaussian distributed with a standard devi ation of 
Q-yx ~ .18 that matches the observational results of Man tz et aP 
( l2010ah . which report a scatter of cr,,f = 0.185. 

In the middle panel of Figure |2] we cor npare the Yx - M50 
relation of our sample to observational data (iMantz et al.l2010al) . 
The model agrees nicely at the high-mass end, but underpredicts 
the observed scaling at low masses by about 20% (at the 1-cr 
level). This is the same level of deviation from the data as in the 
case of Lx, which is more significant due to the smaller scatter 
in the Yx relation. The differential contribution to the thermal 
energy per logarithmic interval in radius (and hence to the inte- 
grated Compton-y parameter) is given by dF/dlogr oc r-'Pth('")> 
with the thermal gas pressure Pth = ngas^B?^- It peaks at scales 
slightly sr naller than R^m with 1 -cr contributions extending out 
to 37?5oo dBattaglia et al.ll20Toh . Hence, the observational scal- 
ing constrains our model on those large scales, quite comple- 



' Scatter is calculated as o-y, = [K,— (A + S Xj)]-} IN - 1 

where the sum extends over the data points Xj, 7,, and A and B are the 
fit parameters. 



mentary to the X-ray luminosity. The deviations at small masses 
either indicates different assumptions about /gas, the gas temper- 
ature, or different selection effects of either observat ional sample 
that we use for model calibration and comparison. iMantz et al] 
(1201 Oah determine their masses by adopting a constant value 
for /gas, in contrast to our appr oach which adopts the observed 
/gas.soo - M50 relation giv e n bvlSun et al.l (2009). Additionally, 
we adopt the i Mantz etaP (1201 Oal) centrally included t empera- 
ture throughout all our work, while iMantz et al.l (1201 Oah use the 
centrally excised temperature to calculate Yx- This assumption 
also impacts the scatter of the Yx- M relation. In fact, using the 
centrally included temperature, we found a scatter of cr,, j ~ 0. 1 1 
(see Table[3]where our Yx scaling relations are reported), signifi- 
cantly higher than the value of CTyx - 0.052 found bv Mantz et aP 
(l20T0ah . 

In the right panel of Figure |2] we compare the Fsz - 
M500 relation in our model (calculated as in Equation (3) of 
Battaglia et a l.il201lh with t he ob served scaling relation by the 
Planck Collab oration et aP (1201 ih . Their sample contains clus- 
ters up to z ~ 0.45 and has a median of z ~ 0.15; hence, we com- 
pare the data to our relation at the MultiDark snapshot z = 0.16 
(that contains 11419 clusters above our mass cut) which how- 
ever is not used throughout the rest of the work. Our model re- 
produces the data remarkably well, except for the high-mass end 
where our simulations have a weaker constraining power due 
to the smaller box size in comparison to the survey volume of 
Planck. We find a scatter of c^y^- ~ 0.11 which compares well 
with the Planck result of cr,, , ~ 0. 1 . In Table |4] we report our SZ 
scaling relations for different redshifts. 



2.4. X-ray luminosity function 

Studies of the XLF got out of fashion during the last years due 
to the difficulties of using the X-ray luminosity for cosmological 
purposes. The X-ray emissivity scales with the square of the gas 
density, which makes it subject to density variations and clump- 
ing. This implies large scatter that causes a large Malmquist bias 
and underlines the necessity of careful mock surveys that need 
to address all systematics. 
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Table 2. Lboi - M500 scaling relations. 



redshift z 


A 






B 


CTyx 





-21.41 ± 


0.11 


1.50 


±0.01 


0.179 


0.1 


-21.30 ± 


0.12 


1.50 


±0.01 


0.179 


0.2 


-21.50 ± 


0.13 


1.51 


±0.01 


0.178 


0.4 


-21.13± 


0.17 


1.49 


±0.01 


0.178 


0.61 


-21.59 ± 


0.22 


1.53 


±0.01 


0.177 


0.78 


-20.73 ± 


0.29 


1.48 


±0.02 


0.177 


1 


-20.45 ± 


0.42 


1.46 


±0.03 


0.177 



Note. Scaling relations are reported in the form of 
logio (Lboi / E(z) h-l 10^* erg s"') =A + B logio M500 / ^70 Mo). 
The relation scatter CTy^ is also shown. 

Table 3. lx,5{)0 - A/500 scaling relations. 



redshift z 


A 






B 







-9.18 ± 


0.07 


1.61 


±0.01 


0.109 


0.1 


-8.85 ± 


0.07 


1.59 


± 0.01 


0.109 


0.2 


-8.82 ± 


0.08 


1.59 


± 0.01 


0.109 


0.4 


-8.79 ± 


0.10 


1.59 


±0.01 


0.108 


0.61 


-8.65 ± 


0.14 


1.59 


±0.01 


0.109 


0.78 


-8.36 ± 


0.18 


1.57 


± 0.01 


0.109 


1 


-8.28 ± 


0.26 


1.57 


± 0.02 


0.109 



Note. Scaling relations are reported in the form of 
logio iEiz) rx.50o / Mo keV) =A + B logu, (£(z) M500 / /j^' ^o)- 
The relation scatter cr,., is also shown. 



Table 4. Fsz,500 - A^soo scaling relations. 



redshift z 


A 






B 







-27.93 ± 


0.07 


1.60 


±0.01 


0.109 


0.1 


-27.74 ± 


0.07 


1.59 


±0.01 


0.109 


0.16 


-27.76 ± 


0.08 


1.59 


±0.01 


0.109 


0.2 


-27.65 ± 


0.08 


1.59 


±0.01 


0.109 


0.4 


-27.57 ± 


0.10 


1.59 


±0.01 


0.108 


0.61 


-27.50 ± 


0.13 


1.59 


±0.01 


0.109 


0.78 


-27.15 ± 


0.18 


1.58 


±0.01 


0.109 


1 


-27.01 ± 


0.26 


1.58 


±0.02 


0.109 



Note. Scaling relations are reported in the form of 
logio {EizT^I^ Ksz.500 / Mpc^) = A+B logio iE{z) M500 / h-'o ^o)- 
The relation scatter cr,., is also shown. 



Nevertheless, it provides a complementary check for our 
model. To this end, we use the ROSA T brightest cluster 
sample (BCS) XLF (lEbeling et al.l Il997h . which is in good 
agreement with results from the ROSAT ESO Flux-Limited 
X-ray (REFLEX; Bohringer et al. 2002) and HIFLUGCS 
jReiprich & Bohringei1l2002 ). Note that the XLF is fully deter- 
mined by the mass function and the Lx - M500 relation after tak- 
ing into account the observational biases. This means that apply- 
ing the Malmquist and Eddington bias corrected Lx - M500 rela- 
tionby Mantz et al. (2010a) directly to the MultiDark mass func- 
tion and accounting for the observational scatter in Lx - M500 
should yield an unbiased XLF. We show the resulting bolomet- 
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0.0 0.5 1.0 1.5 

log,„ [hi lO-"* erg s '] 

Fig. 3. Bolometric and soft-band (0.1 - 2.4 keV) XLFs. Shown 
are the soft-band data points, the soft- band and bol ometric 
Schechter fits of the BCS sample of EbeUn g et all (|1997|) . which 
has a median of z ~ 0.08. While the so ft-band XLF, which was 
obtained applying the lMantz et al.l(l2010ai) Lx - M500 relation to 
the MultiDark z = 0.1 snapshot, compares well with the BCS 
data points, it deviates from the correspondi ng Schechte r fit. 
We also show the bolometric XLF of Mantz et alJ (I2010al) and 
bolometric XLF of our model at z - 0.1. The XLFs are calcu- 
lated in equally log-spaced mass bins; the error bars represent 
the Poissonian errors. Note that we limit the comparison to the 
luminosity range covered by our sample, where we cut the low- 
est part because in that range the XLF rapidly drops due to the 
imposed mass cut. 



ric and soft-band (0.1 - 2.4 keV) XLF in Figure[3]and compare 
those to the coiTesponding BCS XLFs and to our model predic- 
tions. Note that there is only the Schechter fit avail able for the 
BCS bo lometric XLF. While the soft-band XLF by Mantz et alj 
(l20T0ah agrees well with the BCS data points, it deviates from 
the corresponding Schechter fit at low luminosities. This is also 
true in the bolometric band, where the XLFs of jMantzeTal] 
(l20T0ah and our model agree well, but deviate from the BCS 
Schechter fit at low luminosities. This may be an artifact due to 
the use of Schechter fit instead of the data points or may point 
to incompleteness of the BCS sample. Note that the Poissonian 
errors of the XLF obtained from the MultiDark simulation are a 
lower limit as we are neglecting the uncertainty due to cosmic 
variance. Studies of the XLF will become again an important 
topi c with the upcoming la unch of the eROSITA satellite (see, 
e.g.. lCappelluti et ani201lh and further studies in this direction 
are desirable. For these reasons, we do not show XLF predictions 
at other redshifts, leaving this for a future study. 

Summarizing, our phenomenological approach provides vi- 
able gas densities that reproduce the observed scaling relations 
of the ICM as well as the XLF. Thus, we will apply it in the 
following to model the CR population in galaxy clusters and to 
predict the radio and gamma-ray emission of our sample. 

2.5. Cosmic ray modeling 

We assume a power-law for the CR proton distribution, 
f{R,p)dp = C{R)p ''dp, which is the effective one-dimensional 
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momentum distribution (assuming isotropy in momentum 
space). We are interested in calculating the radio synchrotron 
(and gamma-ray) emission resulting from secondaries of 
hadronic CR interactions with protons of the ICM. To start, we 
provide the synchrotron emissivity /„, at frequency v and per 
steradian, of a steady-state electron population where radiative 
cooling balances injection from hadronic interactions (adapted 
from Pfrommer et al. 2008 and EnBlin et al. 2 01 L) . 



jy = AyC{R)p^,,(R)- 



(a-2)/4 



■ 1 ( ' ^5) 

fB(") + eCMB \ ffic / 

where the abbreviations Ay and eg^ are defined in Appendix IB] 
fCMB is the energy density of the cosmic microwave background 
(CMB), and eg - l{%n) denotes the magnetic energy density. 
We assume a scaling of the magnetic field with gas density that 
is given by 



B{R) = Bo 



Pgas(^) 

Pgas.O 



(6) 



where Bq is the central magnetic field and ob describes the de- 
clining rate of the magnetic field strength toward the cluster out- 
skirts. Such a parametrization is suggested by cosmological sim- 
ulations (Dubois & Tevssier 2008) as well as Faraday rotation 
measurements (B onafede et al., ,2010l: iKuchar & EnBlinI 1201 iL 
and references therein). The radio surface brightness Sy(R_i_) (in 
the small-angle approximation) and luminosity Ly, at a given fre- 
quency V, are given by 



J = 2 



R 



zdR, 



Ly = 4n dVjyiR) 



(7) 
(8) 



The flux is given by Fy - LyjiAnD^) where D is the luminosity 
distance to the object. Note that we do not convolve S y with the 
instrumental point spread function unless specified. 

The spatial CR distribution within galaxy clusters is gov- 
erned by an interplay of advection, CR streaming and diffusion. 
The advection of CRs by turbulent gas motions is dominated by 
the largest eddy turnover time Tm ~ Ltu/^tu- Iau denotes the tur- 
bulent injection scale (typically of order the core radius) and i/tu 
is the associated turbulent velocity that gets close to the sound 
speed Vs for transsonic turbulence after a cluster merger and re- 
laxes to small velocities afterwards. In contrast, the crossing time 
of streaming CRs over Ltu is ~ Xb ^tu/vst with the streaming 
velocity given by v^^ ~ Vs and Xb ^ ^ parametrizes the mag- 
netic bending scale. Magnetic bottlenecks for the macroscopic, 
diffusive CR transport, are critical in lowering the microscopic 
streaming velocity of CR by some finite factor. Therefore, we 
can define a turbulent propagation parameter 



rtu = 



Tst Xb Vtu 



Ttu 



Vst 



(9) 



that indicates the relative importance of advection versus CR 
streaming as the dominant CR transport mechanism. After a 
merger, turbulent advective transport dominates yielding ytu s> 
1, which results in centrally enhanced CR profiles. In contrast, in 
a relaxed cluster, CR streaming should be the dominant transport 
mechanism implying ytu ~ 1 and producing flat CR profiles (fo r 
a detailed discussion of these processes, see lEnBlin et al.ll2bl ih . 

What we propose here is to keep the spectral shape of the 
CR distribution function as determi ned from cosmological hy - 
drodynamical simulation of clusters (iPinzke & Pfrommerl2010l) . 



which however did not account for CR streaming. Hence the spa- 
tial CR distribution has to be modified to include the effects of 
CR streamin g. For this reason we adopt the analytical result from 
En fiUn et al.l (|201 li) . This yields a model that includes the nec- 
essary CR transport physics and is able to predict the radio and 
gamma-ray emission while reproducing the observed RH prop- 
erties. Note that this approach is not fully self-consistent and 
points to the necessity of future hydrodynamical simulations to 
include the effect of CR streaming and diffiision on the CR spec- 
trum. 

To construct su ch a model, we hav e to generalize the ap- 
proach proposed by lEnBlin etanJMT ) in order to account for 
GNFW gas profiles (Section [2.2l l. We also have to include the 
cluster m ass-scaling of th e CR normalization obtained from sim- 
ulations dPinzke & Pfrommerii2010i) . While details are given in 
Appendix |A] we show below the main steps. When turbulent 
advection completely dominates the CR t ransport, the CR nor- 
malization can be written as (.Enfilin et al.ll201 li) 



C.dv(/?) = Co ^ ' ^QMRf" 

\ ^th,0 / 



(10) 



where jScr = (a -h 2)/3, y = 5/3. We introduced the advective 
CR profi le t](R) - (P(^)/Po) ^^y. Solving the continuity equation 
for CRs. lEnBlin et all (1201 ih derive the CR density profile. 



Pcr(R) = pcR,o?7(^)exp 



(11) 



where /?, - ytu^c and R^ is the characteristic radius, of order the 
core radius, at which the turbulence is supposed to be injected. 
Now, we introduce the semi-analytical mass-dependent normal- 
ization of the CR profile of Pinzke & Pf r ommer. (.20 10.) such that 



7,(7?) = 



Cadv(^) 
Co 



ChybndiR) 



I /He 



(12) 



which effectively redefines Cadv(^) by that of our hybrid model, 
i.e.. 



Chybrid(^) - C(R) 



Pg^R) T(R) 
nip To 



(13) 



H ere, C(R) is the norm alization CR profile of Equation (22) 
of lPinzke & Pfromm ei^ (2010). We additionally account for the 
temperature decline toward the cluster periphery, T(R), given by 
the fit to the universal temperature pro file obtained from cos- 
mological hydrodyn amical simulations dPfrommer et alJ l2007t 
Pinzke & Pfrommer 2010) and deep Chandra X-ray observa- 
tions ( Vikhlinin et al. 2005 ). Eventually, the CR profile in our 
hybrid model is C{R) = Co(pcr(/^)/pcr,o)^™ within R^ of 
Equation ( IA.2l i. with pcR defined by Equation ( fTTT i where Chybrfd 
enters through our redefinition of rf, and C(R) - C{R±) for 
R > R+ and R < R , respectively. 

The last step is to generalize the case of one CR population 
with a single spe ctral index a to include t he spectral curvature 
as suggested bv lPinzke & Pfrommed (l2010l) . They model the CR 
spectrum with three different power-law CR populations with 
spectral indices of a, = (2.15,2.3,2.55). Our formalism can be 
easily extended to account for multiple CR populations by ex- 
tending the terms with a single a to su ms over the three spectral 
indices ('see lPinzke & Pfrommeil20lol) . This extension has to be 
applied to Ay, to the factor (€b(R) / esj''"'^^^'^ of Equation and 
to Equation ( fTOb . While the modifications are straight forward in 
the first two cases (and are adopted for our hybrid model, see 
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Fig. 4. Left panel. We show our hybrid model profiles for the normalization of the CR distribution for the NCCC and CCC cases 
and for different values of ytu- We fix the CR number for the case of ytu = 100 using Equation (36) of lEnBlin et al] (|20TT1) . while 
integrating the cluster volume within Rim, and require CR number conservation during CR streaming. Right panel. We compare 
the hybrid model (adopting ytu - 100) with the semi-analytical ad vection-only case (adopti ng our GNFW gas profiles and the outer 
temp erature decre a se to the simulation-derived model proposed bv lPinzke & Pfrommerl201 0) and with the exact analytical solution 
as in lEnBUn et alJ (120111 but for our GNFW profiles (adopting a - 2.3 and ytu — 100). Here, the CR profiles are normalized at 

Co.l = C(0.1iJ200). 



AppendixlB), introducing a sum over a, in Equation (fTOl l would 
make impossible to solve analytically for rjiR) in Equation ([12 
For simplicity, we decided to only use a - 2.3 in this last case 
For the highly turbulent cases, i.e., for ytu = 100 (1000), we re- 
cover the radial shape and normaliza tion of the semi-analytical 
model of lPinzke & Pfrommed(l2010h within 1% (0.1%). 

Summarizing, our hybrid model for the CR distribution func- 
tion, has the following properties: it accounts for (i) the X-ray- 
inferred (CCC and NCCC) gas profiles and cluster-mass scal- 
ing of the gas fraction, in addition to the universal tempera- 
ture drop in the outskirts of clusters, (ii) a cluster-mass depen- 
dent CR normalization and universal CR spectrum as derived 
from cosmological hydrodynamical simulations, (iii) an effec- 
tive parametrization of active CR transport processes, including 
CR streaming and diffusion, which allows us to explore different 
turbulent states of the clusters in our MultiDark sample. 

In the left two panels of Figure H) we show our hybrid CR 
normalization for the NCCC and CCC cases and for different 
values of ytu- As expected, when CR streaming is the dominant 
CR transport mechanism, i.e., for negligible advective turbulent 
transport or equivalently, ytu ~ 1, the spatial CR profiles are 
flattened irrespective of the cluster state. While turbulence in 
NCCCs could be injected by a merging (sub-)cluster, in the case 
of CCCs, the interaction of the AGN jet or radio lobe with the 
ambient ICM could be the source of turbulence. 

In the right panel of Figure |4] we compare our hybrid model 
profile with the semi-analytical advection-only case (adopting 
our GNFW gas profiles and the outer temperature decrease to 
the model proposed by Pinzke & Pfrommer 2010) and with the 
exact analytical solution as in .Enfilin et al., (.201 li) . but for our 
GNFW profiles (see Appendix lAl for details). The profiles are 
normalized at 0.17? 200- In that case of dominant advective CR 
transport, our hybrid model compares nicely with the semi- 
analytical model derived fr om cosmological cluster simulations 
(iPinzke & Pfrommeill20Tol) . The main differences between our 
hybrid model (and the semi-analytical model) on the one side 



* We checked that the choice of a in Equation dlOt has only a minor 
effect on the results. Varying a within 2.15-2.55 leaves the radial shape 
and the normalization unchanged within 0.5%. 



and the analytical solution on the other side is the inclusion of the 
simulation-based "reference" profile C for the advection-only 
case and the universally observed temperature drop towards the 
outskirts of clusters. Note that the profiles in our hybrid model 
are generally more centrally peaked in comparison to the analyt- 
ical G NFW ca se, which is due to the enhanced radiative cooling 
in the [Pinzke & Pfrommer (2010) simulations that did not ac- 
count for AGN feedback. Thanks to the flexible parametrization 
in our model, this can be easily counteracted by changing y^ and 
Ob, however, at the expense that these parameters are now degen- 
erate with our assumptions on the CR profile in the advection- 
dominated regime and other possible effects that we are not con- 
sidering such as cluster asphericity (see also next Section[3]l. 

3. Radio Surface Brightness IVIodeling 

In this Section we apply our model to reproduce the emission 
characteristics of four well-observed RHs. Hence for the pur- 
pose of this section, we adopt the measured gas and temperature 
profiles derived from X-ray observations of each cluster. Our 
hybrid model includes an overall normalization gcR of the CR 
distribution function and the hadronic ally-induced non-thermal 
emission (see Appendix iBt. Note that this parameter can be in- 
terpreted as a functional that depends on the maximum CR ac- 
celeration efficiency, g(^p,max), (Pin zke & Pfromme3l2010h but 
only for ytu ^ 100. We will additionally study the CR-to-thermal 
pressure Xcr - PcR/Pth, where the CR pressure is given by 



PcR = 



gCRCmpC^ ^ 



i=I 



'1/(1+9-) 



Oi — I 3 — Ui 



(14) 



Here, c is the speed of light, q - 0.8 is the low-momentum cutoff 
of the CR distribution, and the normalization factors of the indi- 
vidual CR populations are given by A, = (0.767,0.143,0.0975) 
CPinzke & Pfrommerii2010i, see also AppendixlBll. 

3. 1 . Modelling individual radio halos 

For our study, we choose the giant ra dio halos of Coma 
([Peiss et all Il997h and Abell 2163 dFeretti et all I2OOII: 
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iMurgia et aP l2009h . both in merging NCCC s, and the ra- 
dio mini-halos of Perseus (Pedlar et al. 1990') and Ophiuchus 
dOovoni et al. 2009; Murgia et al. 2009), both in relaxed CCCs. 
The radio emission of these clusters is representative of a wide 
class of RHs. Additionally, Perseus, Ophiuchus and Coma are 
among the most promising clusters for gamma -ray observations 
(Pinzke & Pfromm er 2010; Pinzke et al.l l201 ih . In Table B we 
summarize the main characteristics of these halos and detail the 
main parameters adopted for the modeling and the correspond- 
ing results. In Figure |5] we show the surface brightness and 
CR-to-thermal pressure profile of each cluster. We model these 
clusters at 1.4 GHz and within 7?2oo- 

3.1.1. The Coma radio halo 

The giant radio halo in Coma has a morphology remarkably 
similar to the extended X-ray thermal bremsstrahlung emis- 
sion, although the radio emiss ion dechne s more slowly to - 
wards the cluster outskirts (Brie l et a n [T99llDeissetan[T997l) . 
Note that the morphology is non-spherical, showing an elonga- 
tion in the East- West direction. The full - width half maxim um 
(FWHM) of the radio beam is 0.156 deg dPeiss et al.lll997h . al- 
most t wo ord ers of magnitude larger than the X-ray observa- 
tion of iBriel et al. ( 1992)0 Thus, we apply a Gaussian smooth- 
ing to our theoretical surface brightness of Equation O with 

0"smoothing = /^W//Mi-adio/2.355. 

We investigate different values for q-b - 
(0.3,0.4,0.5,0.6,0.7) and for ytu ranging from 1 to 100. 
Fir st, we fix the CR nu mber for ytu = 100 using Equation (36) 
of lEnBhn et all (1201 ll) while integrating the cluster volume 
within /?200- Then, we require CR number conservation during 
CR streaming (for CR energies E > GeV where Coulomb 
cooling is negligible for CR protons), i.e., when lowering the 
values of ytu- Fixing the central magnetic field Bq = 5 jiG 
dBonafede et al.ll2010h . we use ^cr as normalization factor to 
match the radio observations. The best match to the data is 
obtained for ytu — 1 and q-b - 0.6. Values as high as ytu ~ 4 
can be accommodated, however, at the expense of a shallower 
decline of the magnetic field profile (smaller as) as a function 
of cluster-centric radius. We recover the radio surface brightness 
shape and the total radio luminosity within about 20%. 

The gamma-ray flux (see Appendix [Cl within /?2oo for the 
ytu = 1 case and for energies above 100 MeV (100 GeV) is 
Fy - 4.1 X 10"^ (1.5 X 10"'^) cm"^ s~'. Hence, it is in tens ion 
with the limit recently set by Fermi-hAI (IZimmer et alJ2012l) of 
Fy.uL(> 100 MeV) =s 2.5 x lO-** cm^^ s'K Note that for sHghtly 
higher values of ytu, i.e., a more centrally concentrated CR distri- 
bution, the radio and gamma-ray yield would be increased (as- 
suming CR number conservation). However, in order to match 
the observed radio synchrotron profiles, we have to decrease the 
CR normalization (parametrized by ^cr)- This causes the asso- 
ciated gamma-ray flux also to be reduced to a level that is low 
enough to easily circumvent the gamma-ray constraints. In fact, 
for ytu = 4, we obtain Fy{> 100 MeV) = 9.6 x lO"'" cm^^ s"' 
(and Fy{> 100 GeV) = 3.6 x 10"'^ cm"^ s"'). We caution that 
in principle, CR streaming should cause the CR spectrum to 
steepen. This may then considerably weaken these constraints 
as a result of the convex spectral curvature since the gamma-ray 
emission probes the high-energy tail of the CR distribution that 
is suppressed in this picture in comparison to the lower-energy 

' Note that the apparent displacement of the radio and X-ray peak of 
about 0.05 deg is well within the angular resolution of the radio obser- 
vation and hence negligible for the modeling. 



proto ns that the radio emission is sensitive to (see lAleksic et al.l 
2012, for an extended discussion of this point). 

No te that this analysis confirms the result by lArlen et al.l 
(120 12h who also conclude that the hadronic model for the 
Coma radio halo is an attractive explanation for magnetic 
field estimates inferr ed by Faraday rotation measure studies 
(Bonafede et al. 2010) and is not challenged by Fermi upper lim- 
its on the gamma-ray emission if the full ran ge of CR spectral in- 
dices is considered, i.e., for 2.1 < a < 2.5. iPfrommer & EnBlinI 
('2004a') demonstrate that those spectral indices provide a fair 
description of the radio emission after taking into account the 
negative flux bowl due to the SZ effect that causes the measured 
radio halo spectrum to experience a sharp cutoff' towards high 
frequencies. 

This is not ne cessarily in conflict with results by 
iBrunetti et al.l (1201 2|) who point out problems of the hadronic 
model for the limited parameter space a > 2.6, while adopt- 
ing a slightly shallowe r azim uthally averaged radio halo pro- 
file ( Brown & Rudnick 201 lb than the one used in this analysis 
(iDeiss et alJl997l) . Careful work is needed to understand the ori- 
gin of the observational discrepancy and to cleanly separate the 
radio halo emission from that of the Coma radio relic and other 
sources in the field. Interestinglv.l Brunetti et al.l(l2012 l) need very 
flat CR profiles both for their hadronic and re-acceleration mod- 
els in order to reproduce the Coma surface brightness which cor- 
responds to our finding. 

3.1 .2. The radio halo in Abell 21 63 

The morphology of the giant radio halo in Abell 2163 is also 
closely correlated to the cluster thermal X-ray structure. As 
in Coma, the radio emission declines towards the cluster out- 
skirts ata_dowerratejn comparison to the thermal X-ray emis- 
sion (iFeretti et alJl200ll) . The morphological appearance is non- 
spherical, with an elongation in the Eas t- West directio n. We use 
the surface brightness map provided by Murgia et al. (2009) for 
which the synthesized radio beam can be approximated by a 
circulai- Gaussian with FWHMadio - 62". Again, FWHMradio 
is larger than the resolution of the ROSAT observation and the 
corresponding gas density profile. Converted to physical scale, 
Csmoothing is of the Order of that in Coma because of the larger 
distance of Abell 2163. Hence, we also apply Gaussian smooth- 
ing. 

We follow the same procedure as in Coma, and adopt a cen- 
tral magnetic field strength of Bq = 5 fiG. The only match to the 
data is obtained for ytu = 1 and - 0.3, i.e., the flattest possi- 
ble surface brightness. We recover the emission shape and the to- 
tal luminosity within about 15%. The corresponding gamma-ray 
flux within /Jaoo is /^y(> 100 MeV) = 4.2 x 10""' cm^^ s^', about 
two orders o f magnitude lower than th e upper limit obtained by 
Fermi-LAT (lAckermann et al.ll2010al ). and Fy(> 100 GeV) = 
1.5 X 10-1^ cm-2 s-i. 

3.1.3. The Perseus radio mini-halo 

The diffuse radio emissi on in Perseus is th e best known exam- 
ple of a radio mini-halo dPedlar et al.lll990l) and Perseus itself is 
amon g the be st studied clusters in X-rays (e.g.. lChurazov et aH 
12003 t iFabian e t al. 

2004|20lll)- 

Also the Perseus radio morphol- 
ogy resembles that in the X-rays. We proceed as before, but now 
adopt a higher central magnetic field strength of Z?o = 10 fiG. 
Such a larger Bo is expected in a CCC with its higher central 
gas density, implying a larger adiabatic compression factor of 
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Table 5. Radio-halo and mini-halo characteristics. 



name 


z 


D 


M 


LlA GHz, obs 


model parameters 


LlA GHz, model 


references 






[h-^Q Mpc] 


[/i^o' Mpc] 


10^' [h-^eigs-^ Hz-'] 




10^' [h-l ergs-' Hz-'] 




Coma 


0.023 


101 


2.15 


0.72 


1,0.6 


0.86 


[1,2, 3] 












4, 0.3 


0.90 




A 2163 


0.203 


962 


2.07 


15.36 


1,0.3 


13.43 


[3, 4] 


Perseus 


0.018 


78 


0.15 


4.40 


3, 0.4 


4.80 


[3, 5, 6] 












100, 0.3 


3.97 




Ophiuchucs 


0.028 


121 


0.41 


0.19 


5, 0.7 


0.19 


[3,4] 












100, 0.3 


0.23 





Note. Top two rows correspond to giant radio halos, while the bottom two rows are radio mini-halos. repre sents the approximate diameter of the 
RHs. Li.4 GHz, obs is the total lumin osity at 1 .4 GHz. Mim andi?200 are taken from'Reiprich & B6hrineei^('2002'). We use X-ray-inferred gas densities 
Peas a nd temperatures for Coma jBriel et aLlll992f) . for A 2163 and Ophiuchus (Reiprich & Bohrinaer 2002), and f or Perseus ([Ch urazov et al] 
12003 ^. For all clusters, we also adopt the outer temperature decrease seen in cosmological cluster simulations fe.g.. lPinzke & Pfrommeri i2010h 
and deep Chandra X-ray observations ( Vikhlinin et al. 2005). Note that despite the fact that Ophiuchus is a CCC with a central temperature dip, 
we adopt a constant central temperature profile because the dip is only important for radii smaller than 30 /jy^ kpc fsee lMillion et alJl20T^ and 
therefore not critical for the surface brightness modeling. The model parameters column indicates the best match values (top row of each cluster) 
and other acceptable values for the modeling (bottom row of eac h cluster), and Li 4 ghz, model is the corres pon ding predicted total luminosit y at 
1.4 GHz within Rjm - S ee the main text for det ails. References: [lllDeiss et al.l ( 119971) [21 iBriel et all ( 119921) [31 lReitirich & Bohringeil d2002l) [4] 
iMureia et al.l ( l2009h [5l lPedlar et J] ( il990l) [6l lChurazov eTS] ( l20(m 



the m agnetic field during the condensation of the cool core (see 
lAle ksic et al. 2010, 2012 for a discussion on the Perseus mag- 
netic field). 

The best match to the data is obtained for ytu - 3 and 
- 0.4, however, values as high as ytu ~ 100, and as low 
as ytu = 2, can be easily accommodated. We recover the sur- 
face brightness profile and the total luminosity within 10%. The 
gamma-ray flux within Rioo for the ytu = 3 case and for energies 
above 100 MeV (100 GeV) is Fy = 1.4 x lO"** (5.1 x lO^'^) 
cm-^ s-'. Adopting y^ - 100 and - 0.3, the corresponding 
gamma-ray flux above 100 MeV (100 GeV) is Fy = 4.9 x lO"'' 
(1.8 X 10"'^) cm"^ s^^ Note that the gamma-ray flux above 
100 IVIeV of the central ga laxy NGC 1275 as measured by 
Fermi is 2 X 10 cm -2 s-i ('Abdo et al.l l2009h . well above our 
model predictions due to hadronically produced diffuse gamma- 
ray emission that is expected to glow from the entire cluster. 

We can compare these predictions with the upper limit above 
1 TeV, and for a region w ithin 0.15 deg aroun d the cluster 
center, recently obtained by lAleksic et al.l (l2012h . For ytu - 3 
(ytu = 100), we obtain a flux of Fy{> 1 TeV,< 0.15 deg) = 
7.3 X lO-i"* (5.5 X lO-'"*) cm-2 s"', which is well below the 
upper limit of the IVIAGIC collaboration, Fy_uL(> 1 TeV, < 
0.15 deg) ^ 1.4 X lO"'-' cm"^ s"'. Note also that, in the case 
of ytu — 100, we obtain a maximum CR acceleration efficiency 
m ultiplier of .?(<rD.max) - 0.52 , about half of the value adopted 
by iPinzke & Pfrommei] (l2010l) . Note that adopting g(^p,max) - 
1 results in slightly smaller gamma-ray luminosities in com- 
parison to those p redicted by Pinzke & Pfrommer (2010) and 
iPinzke et alj (1201 li) because we additionally account for the cen- 
tral temperature dip and as well as the decrease towards larger 
radii. 

3.1.4. J\\e Op[niuc[nus radio mini-[nalo 

The Ophiuchus cluster has been widely studied both in radio 
and X-rays in the last few years because of the claimed presence 
of a non-thermal hard X-ray tail ( Eckert et al. 2008; Fujita et aL 
l2008HGovoni et al.ll2009l:lMurgia et al.„2009ciPerez-Torres et al. 



'2009t iNevalainen etai]|2009 f: 'Murgiaet al.l 12010^, 'Million et aL 
20 1()!). It was classified a s a m erging cluster by Watanabe et aL 
(2001), but more recentlv lFuTita et al.l ( 12008 ) did not find any ev- 
idence of merging and, on the contrary, cla ssified it as one of th e 
hottest clusters with a cool-core (see also iMillion et al.ll2010l) . 
The radio mini-halo morphology displays similarities with the 
thermal X-ray emission. For oiir modejmg we use the surface 
brightness profile provided bv lMurgia et all (120091) . 

We proceed as before, adopting a central magnetic field value 
of Bo = 10 juG. The best match to the data is obtained for ytu - 5 
and ub - 0.7. However, values as high as ytu ~ 100, and as low 
as ytu - 2, can be easily accommodated. We recover the sur- 
face brightness profile and the total luminosity within 20%. The 
gamma-ray flux within /?200 for the ytu = 5 case and for energies 
above 100 MeV (100 GeV) is Fy = 1.2 x 10"'" (4.3 x lO"'"*) 
cm"^ s"'. Adopting ytu = 100 and ub - 0.3, the coiTesponding 
gamma-ray flux above 100 MeV (100 GeV) is Fy = 8.3 x 10"" 
(3.1 X 10"'"*) cm"^ s"'. The gamma-ray flux is, in both cases, 
about two orders of magnitude lower than the upper limit ob- 
tained by Fenni-LAI ( Ackermann et al. 2010a). Note also that 
in the case of ytu = 100 we obtain a maximum CR acceleration 
efficiency multiplier of g(^p,max) = 0.014. 

3.2. Discussion and conclusions on the hybrid hadronic 
model 

After these detailed considerations on four individual RHs, we 
are now summarizing the main points: 

1 . The radio mini-halos in Perseus and Ophiuchus can be easily 
reproduced in our hybrid hadronic model with a wide range 
of possible values for ytu. In contrast, the large extension of 
the giant radio halos in Coma and Abell 2163 in comparison 
to the X-ray emission requires lower values of ytu < 4 in 
order to match the data. However, we would expect higher 
values for yt., in rn erging NCCCs that are highly turbulent 
dEnBlin et al.ll201 11) . There are four possibilities that can ex- 
plain this apparent contradiction, (i) Primary electrons accel- 
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Fig. 5. Surface brightness modeling of the RHs in Coma, Abell 2163, Perseus and Ophiuchus. The left panels show the RHs' 
azimuthally averaged surface brightness profiles, while the right panels show the corresponding CR-to-thermal pressure profiles 
XcR{r) where we additionally mark the RHs' radial extension by a vertical bar. Note that different parameter values that yield 
almost the same surface brightness profile, result in very different Xqr profiles. In the case of Abell 2163 and Ophiuchus, we adopt 
a 10% uncertainty range instead of the errors reported by Murgiaet al. (2009) to account for addition al systematic uncertainties, 
e.g., residual point source contamination. We also adopt a 10% uncertainty range for Perseus, for which jPedlar et al.] d 19901) do not 
report any uncertainty intervals. 
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erated in outer shocks, which we do not consider here, could 
have an important contribution to the RH emission, partic- 
ularly in the cluster outskirts (Pfrommer et al. 2008). While 
this does not have to significantly enhance the total luminos- 
ity, it could account for the larger radial extent of the radio 
surface brightness in comparison to the thermal X-ray emis- 
sion, (ii) Merging clusters are not spherically symmetric as 
can be seen in Coma and Abell 2163, requiring inherently 
non-spherical modelling. In order to reproduce the more ex- 
tended radio emission relative to the thermal X-ray emission, 
the non-thermal clumping factor, Cnon-th, needs to be larger 
than its thermal analogue, Cth, in concentric spherical shells, 
where we defined those statistics by 

Cnon-th - (pgasC) / ( VPgas^) , (15) 
Cth = (p^4/(Pgas)'. (16) 

This appears to be the case since the radio surface brightness 
emission is more elongated than its counterpart in thermal X- 
rays and needs to be carefully quantified in future work, (iii) 
CR streaming and diffusive transport may cause a possible 
spectral steepening in the cluster core region that would im- 
ply spatial variations of the CR spectral index throughout the 
cluster and hence a spatially varying radio emission, in par- 
ticular in the clust er outskirts, (iv) Adopting the simulation- 
derived C profile (Pinzke & Pfrommer 2010) for our hybrid 
model may have biased the inner slope of the CR density 
profile to become too steep due to the overcooling problem 
of purely radiative simulations. This produces cluster cores 
that are too dense (in comparison to observations), which 
also should overestimate the rate of adiabatic compression 
that is experienced by the CR population during the forma- 
tion of the cooling core. Hence, the resulting values of ytu 
are then biased low in comparison to a potentially shallower 
slope of the inner CR profile. To quantify the last point, we 
try to reproduce the Coma surface brightness using a model 
without C. We find that values as high as ytu ~ 8 can be 
accommodated. However, ytu - 1 still represents the best 
match to the data, demonstrating that the problem can be 
weakened but not circumvented even in this case of a cored 
CR profile. Nevertheless, we can successfully reproduce the 
morphology of Coma and Abell 2163 RHs, solving the pre- 
vious problems of th e classical hadronic model (see, e.g., 
iPonnert et al.|20iPjl)- 

2. Gamma-ray and radio emission are complementary tools in 
constraining the CR population. While gamma-ray emission 
is a direct way to constrain the CR pressure, it is also critical 
in disentangling between the hadronic and re-acceleration 
model. In contrast, the information derived from the radio 
window is degenerate with the assumptions on the magnetic 
field strengths. Nevertheless, we find that RH measurements 
allow tighter constraints on the CR pressure (except for the 
case of Coma). 

3. The magnetic field values adopted here are perfectly in 
agreement with other observational constraints, solving pre- 
vious tensions of the classical hadronic model (see, e.g., 
[feltema & Profumo 2011). Indeed, different values of Bq 
could be adopted without violating other observational con- 
straints. The exception to this is again Coma, for which a 
higher Bo value wo uld be in contradiction with Faraday ro- 
tation measure data (Bonafede et alJl2010l) . leaving room for 
central magnetic field stre ngths around 1-5 fiG that are in 
agreement with the data dArlen et al.ll2012l) . Only a value 
much lower would imply a level of the gamma-ray emission 



that c hallenges the most recent Fermi limit (IZimmer et al.l 
I2012h . 

4. Considering the case of advection-dominated CR transport 
(ytu ^ 100), which only allows for a good match to the 
data of Perseus and Ophiuchus, the gcR parameter can be in- 
ter preted as the maxini u m CR acceleration efficiency used 
in iPinzke & Pfromm"erl (12010 ). If the cluster CR popula- 
tion is mainly accelerated in cosmological structure for- 
mation shocks, then this value should depend on the mass 
accretion history and should be approximately universal, 
i.e., similar for all clusters. We find ^cR.Perseus — 

0.52 and 

^CR.Ophiuchus = 0.014, bccausc we fixed Bo = 10 ;uG in both 
cases and used gcR as normalization. This discrepancy can 
be resolved by increasing/lowering the central magnetic field 
in Perseus/Ophiuchus to Bo.Pei-seus ~ 20 fiG and Bo.Ophiuchus ~ 
1 juG. We note, however, that without the guidance of cos- 
mological cluster simulations that include CR streaming, the 
data does not yet constrain ytu- 

5. The small cluster sample analysed here was only meant to 
serve as a proof of concept and to demonstrate the viabil- 
ity of matching observed RH data with our hybrid hadronic 
model. This is in especially true because of the compara- 
bly large parameter space that the physics modeled here de- 
mands and which encompasses our parametrization of the 
magnetic field (Bq, as), ytu, and the CR acceleration effi- 
ciency function. Particularly, a large region of the (ytu, ^b) 
parameter space is available for the analyzed radio mini- 
halos (note also that the parameters ytu and are degen- 
erate). Future work on all known RHs may provide better 
insight, but is beyond the scope of this work. 

Summarizing, we find that the hybrid hadronic model that in- 
cludes microscopic CR transport successfully reproduces the 
main morphological characteristics of both giant and mini-halos 
without violating gamma-ray constraints. As we will see in the 
next section, it can also reproduce the radio-to-X-ray scaling re- 
lation, solving the previous problems of the classical hadronic 
model, and the radio-to-SZ scaling relation. 

4. Radio Scaling Relations 

As introduced in Section [1] there exist an apparent bimodal- 
ity between the radio and X-ray cluster emission. Clusters with 
a given X-ray luminosity can either host RHs or show an ab- 
sence of diffuse ra dio emission (see, e.g., IBrunetti et al.l 120091: 
lEnBUn et al.l 1201 ll) . More recently, a study of the radio-to-SZ 
scaling relation revealed the absence of a strong bimodality di- 
viding the cluster population into radio-loud and radio-quiet 
clusters (Basu 2012). Since Fsz correlates more tightly with 
cluster mass than Lx, this may indicate that the larger scatter 
of Lx correlates with the scatter of the radio luminosity in such 
a way that it produces a bimodality; but as a result of a second 
(hidden parameter) rather than the fundamental cluster parame- 
ter - its mass. In this section, we investigate these two scaling 
relations in the framework of our extended hadronic scenario. 

4. 1 . Exploring the parameter space of scaling relations 

In Figure |6] we show the general scaling relations of our hy- 
brid CR model of Section 12.51 applied to the MultiDark sam- 
ple. We show how both the radio-to-X-ray and the radio-to-SZ 
scaling relations differ upon varying the parameters ytu, Bo, cb 
and redshift. We fix the CR-normalization parameter gcR to 0.5 
in all cases, ensuring an average CR-to-thermal pressure of 2% 



12 



Zandanel, Pfrommer & Prada: On the Physics of Radio Halos in Galaxy Clusters 




Fig. 6. Radio-to-X-ray and radio-to-SZ scaling relations as predicted by our hybrid CR model. In the left panel, we show how the 
LiA GHz-i'X.boi relation varies upon changing different parameters. In the right panel, we show the same, but for the L1.4 ghz-J'sz,500 
relation. Note that in each plot there are two separated populations for each model realization, shown with the same color. The upper 
sets of points correspond to the CCC population while the lower sets correspond to NCCCs. The plot labels indicate those parameters 
which are kept fixed. We fix the gcR-normalization parameter to 0.5 for all cases. See main text for details. 



(0.05%) within R^qo (R500/2). Here, the radio luminosity is cal- 
culated at 1.4 GHz within /?5oo0 In our CR model, we fix the 
CR n umber for ytu = 100 using Equation (36) of Enfil in et al.l 
(1201 1). integrating up to R500. To compute the radio luminosity 
for different values of ytu, we employ CR number conservation 
(for CR energies E > GeV where Coulomb cooling is negligible 
for CR protons). 

In each panel in Figure |6] there are two separated popula- 
tions for each model realization (i.e., for a given set of param- 
eters). Each upper set of points corresponds to the CCC popu- 
lation while the lower set corresponds to NCCCs, respectively. 
In our model, the radio and X-ray emissivities scale with the 
square of the gas density so that L1.4 ghz and Lx.boi are signif- 
icantly higher for CCCs in comparison to NCCCs. In contrast, 
Ysz only depends weakly on the central gas density as explained 
in Section I23] This explains the relative location of the NCCC 
and CCC populations in the L1.4 ghz - Lxm and L1.4 ghz - Ysz 
planes. In particular, CCCs are shifted to the upper right in the 
L\A GHz - i-x.boi plane while they are shifted vertically upward in 
the plane spanned by Li 4 ghz - i'sz- In reality, we expect an (ab 
initio unknown) distribution of these parameters which would 
substantially increase the scatter in the scaling relations and pos- 
sibly lead to a bimodality, depending on correlations among the 
different parameters. 

Most interestingly, the slope of the radio scaling relations 
does not differ when varying parameter values because we do 
not include any cluster mass-dependence in our parametriza- 
tions which is neither required nor constrained by current data. 
Closely inspecting Figure |6] we see that we obtain the largest 
changes in Li 4 ghz for variations in 1 < ytu < 5 and Bq over the 



The mean (median) difference between calculating L,, within R200 
ori?50ois5.3% (5.6%). 



parameter range probed, albeit with a stronger dependence for 
weaker field strengths (as expected from Equation|5]l- 

4.2. Comparison to observations 

After collecting the X-ray luminosity and the SZ flux of known 
RHs, we compare the resulting scaling relations to a phenomeno- 
logical model realization that was chosen to additionally obey 
other observational constraints (e.g., from Faraday rotation stud- 
ies) as well as theoretical considerations on CR transport. 

4.2.1. Observational samples 

In AppendixiDjwe construct a sample of giant radio halos (black) 
and radio mini-hal os (red), a s well as upper hmits on the radio 
emission (Brunetti et al.l 1200 9: Enfilin et a l1 l201lHGovoniet al.1 
2009), and show this in the left panel of Figure I?] The me- 
dian redshift of this sample is z a; 0.18. The corresponding ob- 
servational scaling relation is well fit by logmLi 4 ghz = A + 
B logioLx.boi withA = -37.204+ 1.838 and B = 1.512±0.041, 
and a scatter of o-y^ ^ 0.52 (we do not include upper lim- 
its in the fit; u nits a re as in Figure |7j. We r efer the reader to 
iBrunetti et al] (|2009|) and lEnlSlin et all (1201 ll) for an extensive 
discussion on this topic. We emphasize that in contrast to gi- 
ant radio halos, mini-halos span a wider range in radio lumi- 
nosity (as also pointed out by Murgia et al. 2009). The Perseus 
mini-halo (highest radio mini-halo luminosity in the left panel 
of Figure |7]i, e.g., has a radio luminosity that is almost an order 
of magnitude higher than in giant radio halos at the same X-ray 
luminosity. In contrast, the Ophiuchus mini-halo (lowest radio 
mini-halo luminosity in the left panel of Figure|2ll, which is rep- 
resentative of a few other similar examples recently detected in 
CCCs (such as A2029 and A1835), has a radio luminosity which 
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Fig. 7. Radio-to-X-ray and radio-to-SZ scaling relations in our hybrid CR model (see main text for the details of the chosen parame- 
ters) compared with observations. Left. Li 4 ghz - ^'X.boi relation in comparison to the observational sample take n from the li terature 
and detailed in AppendixlP] Right. Li 4 ghz - J'sz relation in comparison to the observational sample C07 from lBasul (l2012h . 



is much lower than giant radio halos in merging clusters and is 
even below the upper limits. 

We caution that the determination of the slope of the observa- 
tional Li 4-Lx,boi relation is not very robust because of the small 
sample size of RHs, selection biases of extended low-surface 
brightness objects, and systematic uncertainties in the measure- 
ments of Li 4 and Lx,boi- The recently detected low-luminosity 
mini-halos exemplify those uncertainties and the large intrinsic 
scatter of this relation. On the other hand. X-ray luminosities for 
the same object as derived by, e^., ROSAT and Chandra can 
easily differ by a factor of a few F^In the left pan el of Figure |7] 
we additionally show the model of iKushnir et al.l (|2009) with a 
slope of a; 1 .2, arbitrarily normalized for visual purposes, from 
their simple analytical hadronic model. 

In order to compare our model to the obser ved 1.4 GHz 
radio-to-SZ scaling relation, we use the result by 'Basu' ('2012') 
that is based on the radio halo sample from Cassano et al. (2007) 
(hereafter C07), which has a median reds hift of z ~ .18 and ex- 
cludes radio mini-halos. For this sample, |Basu| (|20 12i) compute s 
Ysz within the radio-halo radii given by Cassano et al.l (l2007h . 
These radii have a median of about 0.5 /i^q Mpc which com- 
pares favourably with the median of R^oo x 0.4 /i^q M pc in 
our M ultiDark z - 0.2 snapshot. For the C07 sample, iBasul 
obtains a scaling relation in the form of logjg Li 4 ghz - 
A + B logio Ysz with A = 29.7 + 0.8 and B = 1.17 + 0.18, and 
a scatter of c^y^- ~ 0.28 (units are as in Figure |7]l. The same 
comments regarding the small sample size of RHs, selection 
biases, and systematic uncertainties in the luminosity measure- 
ments also apply here. 

4.2.2. Model realization 

In order to compare our model prediction with observations, we 
select a particular realization of our model. To this end, we use 
the MultiDark cluster sample at z = 0.2, which compares well 
with the redshift of the observational samples; see above and 



" For example, the bolometric X-ray lu minosity of A 2163 a s mea- 
sured by ROSAT is 8.65 x lO"*^ h'^ erg s"' dBrunetti et alll2009l) while 
the Chandra measurement is 4.93 x 10*' h:jg erg s"' dCavaenolo et alj 
[2009; ACCEPT: Archive of Chandra Cluster Entropy Profile Tables; 
|http://www.pa. msu.edu/as tro/MC2/accept/^ . 



AppendixlP] We divide our cluster sample randomly into radio- 
quiet and radio-loud clusters, assuming a ratio of 10% of the 
latter We use the turbulent propagation parameter ytu to sepa- 
rate both populations. In radio-quiet clusters, we assign ytu = 1, 
and in radio-loud clusters, we adopt randomly and uniformly ytu 
values in the intervals [40, 80] and [1,5] for NCCCs and CCCs, 
respectively. We note that those values for ytu would produce 
slightly too concentrated surface brightness profiles for giant ra- 
dio halos. As laid out in detail in Section [J!2l this may be an ar- 
tifact caused by a combination of four effects, namely additional 
primary radio emission in the cluster periphery, an enhanced 
non-thermal clumping factor, spatial variations of the CR spec- 
tral index, and our choice of simulation-derived CR profiles that 
may be too concentrated as a result of the overcooling in radia- 
tive simulations without AGN feedback. Moreover, we demon- 
strate in Figure |6]that the normalization of the RH luminosity is 
highly degenerate with our other parameters so that a different 
combination of those could mimic a similar normalization for 
smaller values of ytu. 

Our modeling of magnetic fields is inspired by Faraday ro- 
tation studies that point to higher fiel d values in the core re- 
gion of CCCs comp ared to NCCCs dKuchar & EnBlinI 120111; 
Bonafede et al.ll20Tol) . presumably due to the higher adiabatic 
compression factor during the formation of the cooling core. 
Hence, for radio-quiet clusters, we adopt randomly and uni- 
formly distributed values of the central magnetic field Bq in the 
intervals [2.5,5.5] fiG and [5, 10] fiG for NCCCs and CCCs, 
respectively. To account for the potential turbulent dynamo in 
radio-loud objects (characterized by a higher turbulent transport 
parameter in our model), we slightly increase Bq in those ob- 
jects and chose Bq intervals of [4.5,7.5] fiG and [7.5, 12.5] fiG 
for NCCCs and CCCs, respectively. 

We fix Q-B = 0.5 and gcR - 0.5 for all clusters. We note that 
our parameter choices are mostly phenomenologically driven 
with the aim to reproduce observations. The parameter study 
in Figure |6] exemplifies considerable degeneracies so that dif- 
ferent combinations of parameters can potentially result in very 
similar distributions. We emphasize the need of more detailed 
observations of radio halos and in particular of multi-frequency 
correlation studies to constrain the interplay of some of these 
parameters. 
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In Figure |7] we show our model in comparison to the ob- 
served radio-to-X-ray and radio-to-SZ scaling relations. The 
normalization of our model can be arbitrarily varied by chang- 
ing ^CR as long as the result ing Xcr respects the current observa- 
tional constraints (see, e.g. jAleksic et al.ll2012l) and remains be- 
low a few percents. As explained above, our choice of gcR = 0.5 
ensures an average CR-to-thermal pressure of 2% within Rsm, 
allowing for plenty of available parameter space to match the 
observational constraints. Our model is sufficiently flexible to 
either mimic a cluster radio bimodality or not, depending on the 
parameters adopted for the populations of radio-loud and radio- 
quiet objects. 

However, with the given model realization as in Figure |7] 
the separation of the radio-loud and radio-quiet populations is 
substantially larger in the Li 4 ghz - ^x.boi plane than in the 
LiA GHz - Ysz plane, which exhibits almost a continuum distri- 
bution from radio-loud CCCs to the radio-quiet NCCCs. This is 
mainly because the bolometric X-ray emissivity scales with p^^^, 
while Ysz ^ Pgas (which is only strictly valid for an isothermal 
gas distribution). As a consequence, the highly peaked gas pro- 
files of CCCs have less impact on Fsz than on Lx.boi- This is one 
plausible explanation for the observed discrepancy of the pres- 
ence of a bimodality in Li 4 ghz - i-x.boi and the apparent absence 
of it in L1.4 GHz - Ysz- 

The slope of our model depends on the different parame- 
ter choices and, particularly, on the relative differences intro- 
duced for the NCCC/CCC and the radio-loud/quiet populations. 
However, we note that our Li 4 ghz - i'X.boi is slig htly shallower 
than t he observed relation, more similar to the iKushnir et al.l 
(2009) prediction, and that the slope of our L] 4 ghz - Ysz com- 
pares well with the observed slope. Owing to the many uncer- 
tainties and lack of robustness both in the observations and mod- 
eling at this stage, we do not attempt to fine-tune our model to 
the observations. 



5. Radio Luminosity Function 

5.1. Comparison with observations at 1.4 GHz 

In Figure [8] we show the RH luminosity function (RLF) at 
1 .4 GHz for a representative realization of our hybrid CR model 
(as in Section nil, and compare it with observational results. As 
for the XLF, the RLF is completely determined by the cluster 
mass function and the radio luminosity-to-mass relation, through 
LiA GHz - Lx,boi or L1.4 GHz - Ysz in combination with our phe- 
nomenological gas model. However, in the radio band there is 
the additional uncertainty of the fraction of radio-loud clusters. 
Thus, in Figure [8j we also show the RLFs obtained by apply- 
ing the Li 4 GHz - Lx,hoi and Li 4 ghz - Ysz relations to simulated 
clusters of the MultiDark snapshot at z = 0.2, using our phe- 
nomenological gas model for Lx.boi and Ksz.soo of each cluster, 
respectively. Note that this procedure is only applied to halos de- 
fined as radio-loud clusters (which are by definition accounted 
for in the L1.4 ghz scaling relations) and we assume a fraction of 
100%, 25%, 10%, 5% and 1% of radio-loud clusters. As evident 
from Figure [8] this is different from our model for the radio lu- 
minosity, where we also define a fraction of radio-loud clusters. 
However, the radio-quiet population also contributes to the RLF 
with an increasing fraction at low luminosities. This is exempli- 
fied in the top left plot of Figure[8] which shows the contribution 
of radio-quiet and loud populations to the total RLF, assuming a 
fraction of 25% and 1% of radio-loud objects. 

In Appendix |D] we make an attempt to construct an 
RLF from existing X-ray flux-limited radio surveys. There ex- 



ist two such studies, the cluster radio survey done with the 
National Radio Astronomy Observatory (NRAO) Very Large 
Array ( VLA) sky survey (NVSS) at 1.4 GHz of Giovannini et al] 
(Il999l) and the survey with the Giant Metrewave Radi o 
Telescope (GMRT) at 610 MHz bv IVenturi et al] (l2007l l200i . 
For the latter, one can also construct an RLF at 1.4 GHz using 
the corresponding RH follow-up measurements. The fractions 
of radio-loud clusters are about 6%, 18% and 24% for the NVSS 
1.4 GHz, GMRT 610 MHz and GMRT 1.4 GHz samples, respec- 
tively. As explained in Appendix iDl we use the 1.4 GHz NVSS 
RLF (with a median redshift of z x 0.18) as observational ref- 
erence for our comparisons. We conclude that the observational 
determinations of the RLF is not very robust at this stage; the 
very different fractions of radio-loud clusters found in the differ- 
ent studies is one indicator of this. 

Generally, there is fair agreement between the NVSS RLF 
and both our modeled RLF and the RLFs based on observa- 
tional scaling relations, particularly for radio-loud fractions be- 
tween 10% and 1%. In particular, we verified that the cumula- 
tive number of RHs above a certain flux limit of the NVSS sur- 
vey is well matched by our 10%-loud fraction case, which will 
be used in the following section. The RLF obtained from the 
L\A GHz - i'X.bol relation differs from the NVSS RLF at high lu- 
minosities, presumably caused by the large observed scatter. On 
the other side, the RLF obtained from the Li 4 ghz - Ysz relation 
matches the NVSS result better These results need to be con- 
solidated by RLFs corrected for flux-incompleteness and simu- 
lations of larger cosmological volumes that are more complete 
at the high-mass end. Figure [8] shows that it will be difficult to 
discriminate between different scenarios at high radio luminosi- 
ties (or equivalently masses). Indeed, in the bottom right panel of 
Figure [8] we compare our RLF and the RLFs based on observa- 
tional scaling relations for a 10% fraction of radio-loud clusters 
to the NVSS RLF. This suggests that the low-luminosity (low- 
mass) clusters will be the most useful in disentangling between 
different models. This emphasizes the importance of conduct- 
ing homogeneous, well controlled surveys of radio halos with 
the extended VLA at 1.4 GHz and LOFAR at lower frequencies. 
Since the latter has already started to take data, it is extremely 
timely to present RLF predictions for this wavelength regime in 
our model, which we will do next. 

5.2. Low-frequency predictions at 120 MHz 

In Figure |9] we show our model predictions at 120 MHz ob- 
tained with the same representative realization of our model as 
in SectionlH with 10% radio-loud clusters. We show both the dif- 
ferential RLF (top left panel) and the cumulative number density 
(bottom left panel) at different redshifts (corresponding to differ- 
ent MultiDark snapshots of Table 1). We note that our imposed 
mass cutoff of M200 - 10'"^ /i"' Mq, which has been adopted to 
reliably model the cluster gas distribution, translates into a lu- 
minosity cutoff. This causes the differential luminosity function 
to turn over at the low-luminosity end and may artificially flatten 
the slope already at sightly higher luminosities than the luminos- 
ity maximum that indicates our formal incompleteness limit. 

Additionally shown in Figure |9] is the expected LOFAR 
Tier 1 point-source flux limit of F^^ - 0.5 mJy (iRottgering et al.l 
2012) converted to a luminosity limit at a few representative red- 
shifts. This flux limit is clearly an underestimate for nearby RHs, 
which extend over angular scales ~ 1 deg, as e.g., in the case 
of the Coma radio halo. In order to make more reliable predic- 
tions, in the follo wing, we will calcula te the RH flux limit with 
Equation (10) of ICassano et al.l (l2010l) . approximating the halo 
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Fig. 8. RH luminosity function (RLF) at 1.4 GHz. The top left panel shows the RLF of our hybrid CR model (see main text 
for the details of the chosen parameters) for different fractions of radio-loud clusters. Additionally shown is the contribution of 
radio-quiet and radio-loud populations to the total RLF (assuming a fraction of 25% and 1% of radio-loud objects). The top right 
panel shows the RLF obtained by applying the observed Li 4 ghz - i'X.boi relation to the MultiDark clusters at z = 0.2, using our 
phenomenological gas model for Lx.boi of each cluster; again for different percentages of radio-loud clusters. The bottom left panel 
shows the RLF obtained by applying the observed Li 4 ghz - i'sz.soo relation to the MultiDark clusters at z = 0.2, using Fsz,50o of 
our our phenomenological gas model for each cluster. The bottom right panel shows the comparison between the three approaches 
for 10% of radio-loud clusters. In all panels, we show the NVSS survey RLF (Giovannini et al. 1999) with a median redshift of 
z ~ 0. 1 8, corrected for the sample incompleteness and survey sky coverage. Horizontal error bars represent the mass bins while the 
vertical error bars are Poissonian uncertainties. The light gray line marked by the arrow estimates the incompleteness limit owing 
to the adopted low-mass cut (see Section|2]i and the scatter in the halo luminosities. 



extension by Rsqq and requiring the mean flux within the RH 
half-radius to be higher than F^^. This may result in an overesti- 
mation of the flux limit for CCCs, whose radio emission is more 
centrally concentrated than in NCCCs. The median Rsoq of our 
sample is about 0.4 /i^q Mpc at all redshifts. This translates to 
a flux limit of about 30 mJy at z = 0.1, 7 mJy at z = 0.2 and 
0.5 mJy (equal to the point-source value) at z ~ 0.6. 

To derive low-frequency flux functions, we construct an 
analytical model for the evolving RLF. We fit the 120 MHz 
RLF at different redshifts with a second-order polynomial 



of the form logjo dn/dLno mhz = Aq + Ai logio L^o mhz + 
A2 (logio L120 MHz)^- All luminosities are measured in units of 
h^Q&vg s"' Hz"' and (comoving) number densities in units of 

h^Q Mpc""'. We consider only luminosities with log;,) Lua mhz ^ 
29.25 to exclude the turn-over at low luminosities caused by in- 
completeness. To obtain an analytical model for n(L, z), we con- 
strain the evolution of the three free parameters A,- to follow a 
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Fig. 9. RH luminosity function n{L,z) at 120 MHz (top left panel) and cumulative number density of RHs n(> L,z) (bottom left 
panel) at different redshifts z (color coded) for the model realization described in Section |4] with a 10% fraction of radio-loud 
clusters. To obtain an analytical model for n(L,z), we fit the logarithm of the RLF at each z with a second-order polynomial and 
constrain the evolution of the three free parameters to follow a linear function in 1 -H z (see main text for details). The right panels 
show the comparison of the RLF fits (top) and the cumulative number density of the MultiDark samples (bottom) to the constrained 
analytical model. The bottom panels of these two panels show the relative differences An/n = («anaiyticai - nfit)/nfit- Additionally 
shown is the LOFAR Tier 1 point-source flux limit of F^^ — 0.5 mJy (Ro ttgering et al.ll2012h converted to a luminosity limit at a 
given redshift. Horizontal error bars represent the mass bins while the vertical error bars are Poissonian uncertainties. The light gray 
line marked by the arrow estimates the incompleteness limit owing to the adopted low-mass cut (see Section |2]i and the scatter in 
the halo luminosities. 



linear function in 1 -i- z, i.e.. A, = A,,o + A,,i (1 -i- z)[3 In the 
right panels of Figure |9] we compare the RLF fits (top) and the 
cumulative number density in the simulation (bottom) to the an- 
alytical model. In particular the analytics matches the simulation 
well except for high luminosities and high redshifts where small 
number statistics explains the deviations. 

This analytical model describes the number of RHs expected 
in our model per unit luminosity and per unit comoving volume 



The values of these parameters are Ao,o = -436.79, Agj = 109.75, 
A|,o = 29.68, Ai.i = -7.17, A2.0 = -0.51 andAz.i =0.12. 



Vc, i.e., d^N{L,z)/dVcdL. Hence the cumulative number of RHs 
above a given flux limit F is given by the integral 

r- r°° d^N(L,z) dVe 
N(>F)- '.\' -^dzdL, (17) 

Jzi JuF) dVcdL dz 

where L{F) - 47tD{z)^F and D(z) is the luminosity distance 
to an RH at redshift z- The result is shown in the left panel of 
Figure [To] for the model realization described in Section |4] with 
a 10% fraction of radio-loud clusters (black solid line). We limit 
the integral to luminosities logjo L120 mhz ^ 29.25. Our redshift 
integration extends from zi = 0.018, the redshift of the closest 
known RH in Perseus, to zi - 2. As shown in Figure[TO] already 
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Fig. 10. Cumulative number of RHs above a certain flux limit in an all-sky survey at 120 MHz. We show the result of the model real- 
ization described in Section|4]using all clusters and adopting a fraction of 10% of radio-loud clusters (black solid line). Additionally, 
we show the result obtained by using the 10% radio-loud clusters only (black dashed line). We also show the diff'erential contribution 
to the RLF in redshift slices. Note that the number of (detectable) RHs would be dramatically reduced by the presence of a break 
in the model at some low luminosity-scale, or a mass-dependence of the model parameters causing the RH luminosities to decrease 
at low masses. Left. Total number of RHs in the sky. Right. Number of detectable RHs by the LOFAR Tier 1 survey considering its 
sky coverage (half sky) and adopting a realistic flux limit corresponding to diff'erent angular source extensions at different redshifts. 
The point-source flux limit, which applies to distant clusters with z > 0.6, is taken to be F^^ - 0.5 mJy (see main text for details). 



redshifts z > 1 do not significantly contribute to the flux func- 
tion. Additionally, there are large theoretical uncertainties since 
our gas model is not calibrated for these redshifts (and implied 
low-mass range) and very little is observationally known about 
difiiise radio emission on group scales in particularly at these 
redshifts, which motivates our upper redshift limit. 

Figure [10] shows the contribution of different redshift slices 
to the total flux function. We contrast this to the flux function us- 
ing only the subsample of 10% radio-loud clusters (black dashed 
line). This was obtained by constructing the corresponding RLF 
and repeating the steps above in building an analytical model, 
however, discarding luminosities logjg L120 mhz ^ 30.75 in the 
integration. 

In the right panel of Figure [TO] we show the total number 
of RHs that would be detectable by the LOFAR Tier 1 survey, 
where its sky coverage (of about half the entire sky) and the 
signal degradation due to source extensions of close-by RHs is 
taken into account. The latter is calculated with Equation (fTTT i 
and adopting F - F^i^ where Fmin is given by Equation (10) of 
ICassano et al.l (1201 0) as explained above. We clip the flux limit 
for sufficiently distant RHs at the point-source flux limit F^^, 
i.e., for RHs at z > 0.6 we fix Fmin = 0.5 mJy. 

The LOFAR Tier 1 survey at 120 MHz should be able to de- 
tect a total of about 3500 clusters hosting RHs above 0.5 mJy. 
The precise number of detections depends strongly on the un- 
derlying assumptions. There are two main uncertainties in our 
model: the fraction of radio-loud to radio-quiet clusters and the 
corresponding luminosities as well as our assumed RH model- 
ing in low-mass clusters (which are not yet known to host RHs). 
The fraction of radio-loud to radio-quiet clusters is determined 
from a given (degenerate) set of model parameters that include 
7th, Bq, as, gcR, or equivalently Xcr. The fraction of radio-loud 
clusters mostly affects the number of medium-to-high luminos- 
ity RHs (as can be seen from the 1 .4 GHz RLFs in the top left 



panel of Figure |8]l. The total number of RHs is dominated by 
low-luminosity RHs. While this may suggest that the radio-loud 
fraction is of minor importance for the detectable numbers, in- 
stead the opposite is the case. Because of the fixed ffux limit, 
only the most luminous clusters at each redshift are observable 
so that the total number of detectable RHs scales almost linearly 
with the radio-loud fraction. 

We caution that our predicted total number of (detectable) 
RHs depends on the ability to extrapolate the observed and 
modeled scalings down to our adopted mass limit of M200 ~ 
1 .4 X 10'"^ /i^y Mq. If the underlying physics imprinted a charac- 
teristic scale into the CR transport or magnetic field distribution, 
this would manifest itself as a break in the radio luminosity scal- 
ing relations and dramatically reduce (or even increase) the num- 
ber of expected RHs. While this does not interfere with previous 
(high-frequency) measurements at high luminosities, this may 
of critical importance for future, more sensitive (low-frequency) 
measurements of low-redshift clusters that probe the uncertain 
regime of diffuse radio emission in low-mass clusters. 

The reason for this lies in the steep halo mass function which 
ensures that that clusters above a given cutoff (that is either 
physically motivated or observationally realized through a sur- 
vey flux limit) dominate the total number of (detectable) RHs. 
Only if the radio luminosity scaling remains unaltered below the 
survey flux limit, the steepness of the mass function ensures that 
there will be more RHs scattered above the flux limit then below 
it. This is the so-called Eddington bias that causes the inferred 
luminosities based on only the detected sources to end up as an 
overestimate. Hence the presence of any hypothetical break in 
the radio luminosity scaling, which is unconstrained by current 
data, explains the largely varying model predictions in the re- 
cent l iterature, which vary from a few to hundreds observable 
RHs (ICassano et al.ll2010l: ISutter & Ricker 201 j) to thousand s 
of detectable RHs in future surveys (.Enfihn & Rottgeringl2002h . 
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Another relevant issue in such surveys is the identifica- 
tion of RHs and their hosting clusters (see also Cassano et alj 
[2010). RHs constitute a small part of the entire (diffuse as well 
as apparent point-like) radio source population and therefore 
need to be distinguished from the emission produced by other 
sources. A good approach will be to cross-correlate the radio 
maps with high-sensitivity X-ray surveys such as that by the fu- 
ture eROSITA mission, which is expec ted to detect around 1 0^ 
clusters up to redshift z ~ 13 (see, e.g.. lCappelluti et al]|201 ll) . 

Our results show the prospects of the LOFAR survey and 
other future radio instruments in determining the properties over 
a broad range of RH luminosities. In particular, it should permit 
a robust determination of the number of clusters hosting RHs 
at a given luminosity (mass) and to carefully asses complete- 
ness issues. This will be extremely helpful in elucidating the RH 
generation mechanism and in determining the viable parameter 
space for the hadronic model. 



6. Conclusions 

We present a detailed model for the population of RHs in galaxy 
clusters assuming that the radio-emitting electrons are generated 
through hadronic CR interactions with thermal protons of the 
ICM. We improve over previous works by phenomenologically 
exploring the interplay of advective and streaming CR transport 
in modeling the surface brightness emission, RH scaling rela- 
tions and luminos ity functions. W e confirm the theoretical ideas 
suggested by EnBHnetal] (1201 ih and explicitly demonstrate that 
the inclusion of CR streaming is necessary to solve the problems 
of the classical hadronic scenario of radio halos, namely the 
large extents of giant radio halos in comparison to the thermal 
X-ray emission and the bimodality of radio-loud and radio-quiet 
clusters at a fixed X-ray luminosity. 

To this end, we use a complete cosmological sample of 
gala xy clusters selecte d from the MultiDark A^-body simula- 
tion (iPradaet al.ll2"oT2h with redshifts ranging from z = to 1 . 
In a next step, we construct a phenomenological model for the 
cluster-gas distribution. This is characterized by X-ray-inferred 
(CCC a nd NCCC) gas pro files (taken from the REXCESS 
sample, ^Croston et aL 20081) and a cluster-mass-dependent gas 



fraction (Sun et aL r i2009h . We assign such a (cluster mass- 



dependent) gas density profile to each DM halo in our sample 
and sort it into the NCCC/CCC populations according to a dy- 
namical disturbance parameter that is calculated from the DM 
distribution. With this model, we obtain a cosmologically com- 
plete mock catalog of galaxy clusters that matches the observed 
i'X.boi-to-Msoo, ix-to-Msoo, and Fsz-to-Msoo relations, as well as 
the X-ray luminosity function. 

We then construct a new hybrid model for the CR distribu- 
tion in clusters and combine previous results from cosmolog- 
ical simulations and analytical works. We adopt the universal 
spatial and spectral CR distribution found i n hydrodynamic cos- 
molog ical simulations of cluster formation dPinzke & Pfromrnen 
l2010l) . which only account for macroscopic, advective CR trans- 
port. Hence, we supplement this model with the approach of 
[EnBlinet al. (2011) that includes the treatment of microscopic 
CR transport processes. While turbulently-driven CR advection 
can lead to centrally enhanced CR profiles, CR propagation in 
the form of CR streaming and diffusion produces flatter CR pro- 
files, which should be realized for decaying cluster turbulence. 
In our model, we introduce a CR propagation parameter ytu that 
is the ratio of the CR streaming-to-advection time scale. This 
parameter allows us to effectively switch the regimes where ei- 



ther process dominates the CR transport and to explore dififerent 
turbulent states of the clusters. 

This enables us to model the radio surface brightness profiles 
of giant radio halos (as exemplified in Coma and Abell 2163) 
as well as of radio mini-halos (in Perseus and Ophiuchus) at 
1.4 GHz. We demonstrate that our hybrid CR model can suc- 
cessfully reproduce the radio profiles over the entire extents of 
giant and mini radio-halos, solving the problems of the classical 
hadronic model. We recover the total radio luminosity within 
20% for magnetic field configurations favored by Faraday ro- 
tation measurements. Our modeling shows that within the pa- 
rameter ranges that provide acceptable matches to the radio pro- 
files, the resulting gamma-ray emission from the decay of neutral 
pions — an inevitable by-product in hadronic CR interactions — 
is well below observational gamma-ray upper limits provided by 
Fermi and imaging atmospheric Cherenkov telescopes. Future 
and more sensitive gamma-ray observations, with, e.g., the 
planned Cherenkov Telescope Array, may in principle be able 
to disentangle the two competing models for RHs, the hadronic 
and re-acceleration scenarios. 

While we can get excellent matches to the surface brightness 
of radio mini halos for a range of our CR propagation parameter 
of 2 < ytu ^ 100, in the case of giant radio halos, our model 
prefers lower values of 1 < ytu 5 4. However, the morphologi- 
cal disturbances in the X-ray and radio maps strongly suggest a 
highly turbulent state of the ICM in those clusters that host giant 
radio halos. We discuss that a combination of four effects could 
artificially lower the inferred values of ytu^ (i) additional radio 
emission of primary, shock-accelerated electrons in the cluster 
periphery of these mergers, (ii) an enhanced non-thermal clump- 
ing factor in comparison to the thermal clumping factor as sug- 
gested by the larger anisotropy of the radio emission (relative to 
the X-ray emission maps), (iii) a possible spectral steepening in 
the cluster core region due to CR streaming and diffusive trans- 
port that could cause spatial variations of the CR spectral index 
throughout the cluster and hence a spatially varying radio emis- 
sion, in particular in the cluster outskirts, and (iv) our choice of 
simulation-derived CR profiles that may be too concentrated as 
a result of the overcooling in radiative simulations without AGN 
feedback. 

We select a representative realization of our CR model and 
compare it with existing radio scaling relations. Because of CR 
streaming transport and the different gas-density scalings of the 
X-ray luminosity, Lx J pl.^^dV, and the Sunyaev-Zel'dovich 

flux, Y oc J pg.,,kBTdV, our model is able to simultaneously re- 
produce the observed b imodality of radio-loud and radio-quie t 
clusters at the same Lx ('Brunetti et aT]l2009t lEnBlin et al.ll201 ih 
as well as th e unimodal distribution of radio-halo luminosity ver- 
sus Y (lBasull2012l) : thereby suggesting a physical solution to this 
apparent contradiction. We caution however, that some parame- 
ters in our model are degenerate with respect to the resulting ra- 
dio luminosity and radio emission profiles, in particular ytu and 
fffi (our rate of decline of the magnetic field toward the cluster 
outskirts). Multi-frequency data will be needed to better con- 
strain these parameters and to break these degeneracies. 

Assuming a fraction of 10% radio-loud clusters, we demon- 
strate that our mo del matches the NVSS radio halo luminos- 
ity function (RLE. iGiovannini et alJll99"9l) . However, the high- 
luminosity tail of our model RLF is subject to cosmic vari- 
ance because of the comparably small simulation volume of 
Ih^^Gpc^. Interestingly, the RLF derived from the Li.4ghz - 
Lx,bo\ relation differs at high radio luminosities from the NVSS 
RLF; possibly because of selection and incompleteness effects 
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that are not fully taken into account. The comparison between 
different RLFs suggests that the low-luminosity (low-mass) 
regime is the most promising place to differentiate between var- 
ious models. 

It is expected that the next-generation of low-frequency ra- 
dio surveys will probe this regime. Hence, we we make pre- 
diction for the LOFAR cluster survey, in particular, we com- 
pute the 120 MHz RLF and the cumulative RH number den- 
sity. Given our assumption, we would expect the LOFAR Tier 1 
survey at 120 MHz to detect about ~3500 RHs above the flux 
limit of 0.5 mJy. Since the detectable number of RHs is approx- 
imately proportional to the (uncertain) radio-loud fraction, we 
caution that the precise number depends strongly on the under- 
ling assumptions. In particular, we assume that the model pa- 
rameters can be extrapolated down to cluster masses of about 
M200 « 1.4 X 10'"' hjQ Mq without any break in the radio scal- 
ing relation that would indicate additional scales in the physics. 
Most of the RHs in our sample lie at low masses and thus at 
low luminosities that are unconstrained by current observations. 
If, e.g., the magnetic field and/or the CR distribution in clusters 
are not statistically self-similar and would exhibit much reduced 
strength/number density on group scales, our predictions for the 
detectable number of RHs would be dramatically reduced. 

This demonstrates the potential of LOFAR, and other next- 
generation low-sensitivity radio instruments such as APERTIF, 
EVLA and SKA, in determining the RLF properties. In com- 
bination with future X-ray missions like eROSITA, this should 
yield a robust determination of the number of clusters hosting 
RHs at a given luminosity (mass) and thus elucidate the relation 
of the radio emission with the cluster dynamical state. This wiU 
be crucial in determining the parameter values for the hadronic 
model and eventually in elucidating the RH generation mecha- 
nism. 

Summarizing, in this work we have constructed a model for 
the ICM and CR distributions in galaxy clusters that enabled us 
to provide a cosmologically complete multi-frequency mock cat- 
alog for the (non-)thermal cluster emission at different redshifts. 
We will make these catalogs pubUcly and freely available on-line 
through the MultiDark database (www.multidark.org). Those 
contain the quantities pg^s, Lx.boi, Yx, Fsz, ii20 ghz, ghz and 
Ly (among many others). We hope that the community can make 
valuable use of these catalogues in synergy with the future radio. 
X-ray and gamma-ray data. 
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Appendix A: Cosmic Ray IVIodeling Details 

Here, we describe in detail how our hybrid model for the CR 
distribution in galaxy clusters of Se ction | 2.5 1 is constru cted by 
generalizing the analytical results bv lEnBlin et alJ (1201 ih . 



As anticipated in Section 12.51 when advection dominates 
the CR transport, the CR normalization can be expressed as 
in Equation ( fTOl i. However, when CR streaming and diffusion 
dominates, the CR distribution is modified and flattens consid- 
erably. This can be shown analytically by solving the continuity 
equation for CRs and obtaining the CR density profile, pcR, of 
Equation (fTTT i. Assuming P(R)/Po - n^{R)ln(i, i.e., neglecting 
the temperature dependence, and adopting a standard yS-profile 
for the electron density, 



He = no 1 1 + -7 



(A.l) 



lEnBlin et al.l (1201 Ih find that the solution of Equation ( fTTT i is 
physical only for pcR within the radial range R < R < R+ with 



K± - — — Kt 

2y 



1+ A 1- 



2R,r 

3/5,iR, 



(A.2) 



while it is non-stationary outside these radii. In these regions, 
the authors suggest t o set pcr(R) - pcr (R±) for R > R+ and 
R < R , respectivelv. EnBli n et al.l (1201 ih obtain the profile for 
the CR normalization, C(R) = Co(pcR(R)/pcR.of'"', as 



C(R) = Co 



Rl; 



expl— /3cR 



(A.3) 



for R < R < R+, where j6c = 3y8d PcrI'^J, and C{R) = C{R±) 
for R < R and R > R+, respectively. In this way, different CR 
transport cases are parametrized through ytu- A high value of ytu 
characterizes the advection-dominatedc ase, while the CR pro file 
is flat for ytu ~ 1- We refer the reader to lEnBlin et al.l (1201 ih for 
an extensive discussion. 

We want to extend this result to account for (i) our GNFW 
gas profiles of Section 12.21 (ii) the universal temperature 
drop in the cluster outskirts, and (iii) merge it with the uni- 
versal cluster mass-scaling of the CR n ormalization, C, ob- 
tained from hydrodynamical simulations (iPinzke & Pfromrn^ 
,2010) . Therefore, we adopt the hybrid profile, Chybrid = 
C(R)(pg^s(R)/mp)(T(R)/To), of Equation We note that for 
such a choice, there does not any more exist an analytical so- 
lution to Equation (fTTT i as in EnBlin et al. (2011). This results 
in a 5 -order equation and a numerical solution would not be of 
practical use for every cluster in our large MultiDark sample. 
For simplicity, we adopt the formulae above also for our hybrid 
model. I.e., we adopt C(R) = Co(pcr(^)/pcr.{))^™ within R± of 
Equation ( IA.2l i. with pcR defined by Equation (fTTT i where Chybrid 
enters through the advective CR profile, ri{R), of Equation (fTST l. 
and C{R) - C(R±) for R > R+ and R < R , respectively. Note 
that, in our formalism, R^ of Equation ( IA.2b becomes the char- 
acteristic radius of our GNFW gas profile of Section fZ2T i.e., 
Rc - 0.2/?5oo, and ySd = 0.8 (we checked that varying the value 
of ySci between 0.4 and 1.2 has no impact). The relevant factors 
are y^ and the exponential factor of Equation (fTTT i. As we will 
see in the following, the two radii /?± are not critically aff'ected 
by the form of ri{R). Therefore, despite the approximation, this 
approach captures the main CR transport effects. 

Assuming our GNFW gas profile instead of a standard 
;0-profile for the electron density, so adopting P{R)/Po = 
ne,GNFw(^)/no, ther e exist an e xact analytical solution to 
Equation (fTTT i as in lEnBUn et al.1 (f2011t) . In order to evaluate 
the systematic error that we are introducing with the approach 
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where 



1, the effective inelastic cross-section for 



u 



10' 



10° 



10' 



10-' 



NCCC Exact GNFW 
NCCC Approximated GNFW 
CCC Exact GNFW 
CCC Approximated GNFW 



10 




10-" 10-' 



Fig. A.l. Comparison of the exact and approximate solutions for 
C{R) in the case of our GNFW profile. In the latter case, we 
use the formulae by ' EnBlin et"aLr(l201lh but adopt P{R)/Po = 
ne.GNFw(^)/n{) for this comparison only, i.e., we assume here 
an isothermal ICM. We show the GNFW profiles for NCCCs 
and CCCs, as derived in Section 12.21 For each of these profile 
classes, we show three different cases of ytu = 100, 10, and 1 
(top to bottom). We normalize C(R) for the advection dominated 
case (ytu = 100) at R - lO^^R2m and and require CR num- 
ber conservation during CR streaming. Note that the value for 
C(lQ^^R2oo, Ttu = 100) in the CCC case is identical for the exact 
solution and its approximation, while there is a small difference 
of about 9% in the NCCC case. We adopt a = 2.3. 



described above, in Figure lA.ll we compare the exact solu- 
tion for C{R) in the case of our GNFW profile with the ap- 
proxim ate solution where we use the formulae by lEnBlin et al.l 
(1201 ih . and only substitute the electron j6-profile for our GNFW 
profiles (with - 0.2/?50(), jSci - 0.8). In this last case, we 
fix R^ - 10"^/?//?200 to mimic the typical R value of the ex- 
act solution, otherwise an unphysical step feature would appear 
at < 10"^7?/7?2()(). This latter approximation is kept in our hy- 
brid model, which however has no impact on the model surface 
brightness and total luminosity. There is almost no difference 
between the two cases, as shown by Figure lATl This shows that 
the approach presented here to construct our hybrid model can be 
safely followed in order to derive a fully working parametriza- 
tion which captures the main CR transport effects. 



Appendix B: Radio Emission 

The factor Ay in the equation for the r adio synchrotron emis siv- 
ity. E quation (|5]l, is given by (see also lPfrommer et ani2008l and 
lEnfilin etah,2011.) 



16^ ""o-^pinsC^ 

'(tte - 2)0-TeB(,OTp \OTe 



Gev) 



with 



V3^g..3 ,^^VF(^)r(^)F(^) 
32;r m^c^ + 1 



r(^) 



(B.l) 



(B.2) 



proton-proton interactions is cr^ 



pp 



32 (0.96 + e' 



,4.42-2.4» 



Be = ySTrea. ^ 31(v/GHz) fiG, and F is the Gamma-function 
(lAbramowitz & Stegunll965l) . AE,y„j, is given in units of erg, and 
Ay is given in units of erg cm-' g"' sr"' . 

The generalization of the radio emissivity, jy, to account for 
three CR populations, each characterized by different power-law 
spectra an d the inclusion of the norm alization parameter gcR, 
following iPinzke & Pfrommed (1201 Oh . is straight forward. We 
obtain 



jv.hybrid - ^CRC(/?)Pgas(^) 



esiR) + ecMB 



6bW 



(B.3) 



where the sum is over the three CR spectral indexes 
a, - (2.55,2.3,2.15) with the corresponding factors A , = 
(0.767, 0. 143, 0.0975) found bv Pinzke & Pfrommed (l2010l) . 



Appendix C: Gamma-ray Emission 

The gamma-ray flux abov e an energy Ey is given by (see, e.g., 
iPinzke & Pfrommedl2010l) 



FJ> Ey) 



dVAyiR) , 



(C.l) 



where the omnidirectional (i.e., integrated over the An solid an- 
gle) gamma-ray emi ssivity above Ey is Ay(R) = AyC{R)pgas(R)- 
The parameter Ay is dPinzke & Pfromnienl2010h 



Ay - gCRDy,bre^k- 



4m„oc 
3ml 



Si, A,- 



' pp.' 



a/Si \2m„o 
Oi + I a,- - 1 



2Si 



26, 



(C.2) 



where Xj = {1 + [m,ocV(2£y,,)]2'"'), W^Mt = PxM,b) - 
I3x,(a,b) where P denotes the incomplete Beta-function 
('Ab ramowitz & Stegun|[T965h . and 6, = O.UaT^-^ + 0.44. The 
term ZJybreak = Dy{Ey,Ey\,Teik) represents diffusive CR losses 
due to escaping protons from the clu ster at the equivalent pho- 
ton energy for the break fiy break (see Pinzke & Pfronmierll201(]| 
for details). Ay is given in units of cm^ s"' g"'. 



Appendix D: Observational Radio-to-X-ray Scaling 
relation and Luminosity Function 

For comparison with the observed 1.4 GHz radio-to-X-ray 
sca ling relation , we in clude almost all RHs of the sample 
by 'EnBU n et al.l (1201 ll) . We exclude RXCJ 13 14.4-25 15 and 
Z7160 because they lack bolometric X-ray measurements. We 
add to our sample the clusters Ophiuchus, A2029 and A 1835 
(iGovoni et al ] l2009 h. The bolometric X-ray luminosities, Lx.hni , 
of clus ters hosting giant radio halos are taken from Brune tti et al.l 
(120091) . while Lx.boi values of cluster hostin g for mini-halos 
are ta ken from iReiprich & Bohringed (l2002h . iBohringer et al.l 
(fT998l) and ACCEPT. In the cases where measurement un- 
certainties for X-ray and radio luminosities are not reported, 
we adopt a 10% uncertainty. Our final RH sample has a me- 
dian redshift of z ^ 0.18. In the left panel of Figure |7] 
we show the corresponding radio-to-X-ray scaling relation 
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which is well fit by logi()(Li.4 cHz/Ziyo ^rg s ' Hz ') = A + 
B logio(Lx,boi//iTO erg s"'), with A = -37.204 + 1.838, B = 
1.512 + 0.041, and scatter cTyx ~ 0.52. Regarding the clus- 
ters w ith upper lim i ts on the diffuse radio emission in the sam- 
ple of Sin; et al.l (1201 ih . we select only those 8 clusters with 
ACCEPT measurements to obtain a homogeneous data set in 
^x.boi (note that there are a number of clusters with upper limits 
on Li 4 GHz for which there are only soft-band X-ray luminosities 
available). 

In Figure ID. II we make an attempt to construct an RLF 
from existing X-ray flux-limited radio surveys. There exist two 
such studies: the Giovannini et al. ( 1999) survey with NVSS at 
1.4 GHz and the ^Venturi et al... (2007. .2008) survey with GMRT 
at 610 MHz by. We only select RHs, i.e., we do not consider 
radio relics or other difiiise radio emissions of unclear classifi- 
cation. The 1.4 GHz NVSS survey contains 13 RHs out of 205 
analyzed clusters while the 610 MHz GMRT survey contains 
6 RHs out of the 34 obser ved. The sample finally analyzed by 
IVenturi et al.l (l2007l l2008h is composed by 50 clusters and we 
can also build a corresponding RLF at 1.4 GHz using the 12 
present RHs. The fraction of radio-loud clusters is about 6%, 
18% and 24% for the NVSS 1.4 GHz, GMRT 610 MHz and 
GMRT 1.4 GHz sample, respectively. The corresponding me- 
dian redshift is 0.18, 0.26 and 0.25. We calculate the RLF using 
the classical Vmax estimator (see, e.g.. 'F elteiilll976l) correcting 
for the sample incompleteness and survey sky coverage. The 
most problematic aspect in obtaining these RLFs, apart from 
the few available objects, is the calculation of a meaningful flux 
limit. We obtain such a limit by fitting the upper envelope of 
the luminosity-distance distribution of a given sample, as shown 
in the insets of Figure ID. II following the procedure adopted by 
iBroderick et al.l (l2012h . Note that it is particularly ambiguous to 
calculate a meaningful flux limit for the GMRT survey due to 
its poor luminosity-distance distribution. We decide therefore to 
take the 1.4 GHz NVSS RLF as reference in our comparisons 
with observation. However, we want to stress that several issues 
can affect this result, such as the small sample size that impacts 
the resulting flux limit, and the Malmquist and Eddington biases. 
Indeed, the very diff'erent fraction of radio loud clusters obtained 
from diff'erent samples is one clear indicator of the large uncer- 
tainty in the RLF. 
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Fig.D.l. RH luminosity function obtained from existent obser- 
vations. The insets show the sample luminosity-distance distri- 
butions (see main text for details) where the solid line is the fit 
to the upper envelope population (indicated in green) employed 
to calculate the flux limit for the classical Vmax estimator. The 
choice of the upper envelope population is somehow arbitrary, 
particularly in the GMRT cases due to the poor luminosity- 
distance distributions. The horizontal error bars represent the 
mass bins while the vertical error bars are Poissonian uncertain- 
ties. 
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